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Abstract 



A new method of approximation scheme with potential application to 
a general interacting quantum system is presented. The method is non- 
perturbative, self- consistent, systematically improvable and uniformly ap- 
plicable for arbitrary strength of interaction. It thus overcomes the various 
limitations of the exsting methods such as the perturbation theory, the vari- 
ational method, the WKBJ method and other approximation schemes. The 
current method has been successfully applied to a variety of interacting sys- 
tems including the anharmonic/ double- well oscillators ( with quartic-, sextic- 
and octic couplings ) and the scalar field theory with quartic-coupling in the 
symmetric phase. The method yields important insight in to the structure 
and stability of the interacting- vacuum of the theory. The results are in good 
agreement with the exact predictions of supersymmetry where ever applica- 
ble. Possible further applications in the areas of quantum statistics, finite 
temperature field theory, condensed matter physics etc are also discussed. 
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1. Introduction 



It is now generally believed that quantum theory (QT) describes 1 ' 2 the 
observed features of the physical universe at the fundamental level, both 
at the microscopic- and the macro- scales, which span a vast range in dis- 
tance, typically from ~ 1CT 16 cm to light years and time scales, from 
10~ 15 second to the age of the universe ~ 10 10 years. It is truly amazing 
that phenomena as diverse as the colour of the sky, the emission of light by 
heated bodies, stability of matter, the evolution of stars, super conductivity 
and other properties of matter at extreme conditions of temperatures, den- 
sity and pressure etc, all yield to description 2 by the QT. 

However, inspite of the vastness of the range of applicability of QT, the 
exact, analytic predictions of the theory based upon the solution of the under- 
lying fundamental dynamical equations are rather sparse 3 . It is not difficult 
to count and classify the cases where exact analytic solutions are possible. In 
fact, it is currently believed 4 ' 5 that the analytic solvability of the dynamical 
equations of quantum theory is mainly based upon the factorisation property 
of the corresponding quantum-Hamiltonians. In that context, recent applica- 
tions of super-symmetric quantum mechanics (SUSYQM) 6 have resulted in 
considerable extension 6 of the range and class of exactly solvable potentials 
in non-relativistic QT. 

In view of the rather limited range of physical phenomena amenable to 
exact analytic solutions in QT, it has become inevitable to develop various 
approximation methods of solution. This aspect was realized soon after the 
discovery of QT. The history of approximation methods is, therefore, as old 
as the theory itself 7 . Consequently, this has generated a vast amount of lit- 
erature devoted to the various approximation schemes (AS). 

In discussing the merits and limitations of the various AS in QT certain 
basic criteria are generally applied. These criteria may, for example, include 
the following: 

(i) general applicability (ii) simplicity of implementation (iii) achieved 
accuracy (iv) systematic improvability (v) rapidity of convergence and 
adaptability for numerical computation (when the obtained results are not 
in closed formjetc. In spite of the vast amount of accumulated literature 
and the broad range of the existing AS, there is considerable current effort 
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in searching for further improved AS which would satisfy most of the criteria 
listed above. 

The present thesis records original research aimed at achieving a novel 
method of approximation in QT, which is non-perturbative, self-consistent, 
systematically improvable and simple to implement. More importantly, it 
has the potential to be applicable, in principle, to any general and arbitrary 
interacting system including quantum field theory and quantum statistics. 

In order to further motivate and understand the scope of the present 
scheme, it is necessary to review and survey the various existing approxima- 
tion schemes (AS) in QT and discuss their relative merits and limitations 
as against the set of desirable criteria listed above. However, such a task 
is understandably rather arduous due to the vast amount of accumulated 
literature. Moreover, several popular schemes of approximation already find 
mention in texts on QT. We shall, therefore, confine ourselves to the main 
AS's and shall be rather brief. In the following sections a survey along the 
stated lines is presented for the following main categories of the AS in QT : 

(1) Methods based upon Perturbation theory 

(2) The Variation Method of Approximation 

(3) The Hill- determinant and related methods 

(4) Combination of variation method and perturbation techniques 

(5) The Algebraic Method of Operator Expansion 

(6) The JWKB approximation method 

(7) The Method of Approximation employing the Canonical Transforma- 
tion in QT 

(8) The Approximation methods based upon super symmetric quantum me- 
chanics 

(9) The methods based upon the re- summation techniques 
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(10) The Approximation schemes based upon the path-integral formu- 
lation of QT 

(11) The self- consistent schemes of approximation 

(12) Other methods of approximation 

In the next Chapter, we make a brief survey in dealing with each scheme 
of approximation and highlight the main features of each with special em- 
phasis on its merit and limitation. Each method is also applied to the case 
of the quantum an-harmonic oscillator (QAHO) taken as the reference stan- 
dard. This quantum system is chosen for the above purpose since it is also 
the traditional choice as the theoretical laboratory for testing various AS in 
QT 8 ' 9 . Consequently, the QAHO is perhaps among the most widely studied 
quantum system, see, e.g., refs. [8,9]. 



2. Survey of Approximation Schemes in Quan- 
tum Theory 

In the following sub-sections we present a brief survey of the different 
main schemes of approximation in QT as listed in the previous Chapter. 

(1) Methods based upon Perturbation- Theory 

This scheme is among the earliest approximation methods recorded in the 
literature. In the context of classical physics, the method was applied by 
Lord-Rayleigh 7 as early as 1873. In QT the application of the method is 
as old as the theory itself, see, e.g. ref.[7]. Since the salient features of the 
method are well discussed in texts (e.g. [3]) only brief discussion follows. For 
simplicity, we restrict to the case of non-degenerate and discrete spectrum. 

The generic central feature of the perturbation method consists in ad- 
dressing to those problems where the system Hamiltonian can be split in the 
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following fashion 



H = H s + \Hj , 



(1) 



where Hs corresponds to an exactly solvable Hamiltonian whereas the other 
part, AH i can be regarded as a perturbation-interaction term added to the 
dominant Hamiltonian represented by H s , the implicit assumption being 
that the perturbation term is "small" in some-sense ( to be specified be- 
low). In the above equation, the parameter 'A' generically corresponds to 
an expansion-parameter. This may be either naturally present in the inter- 
action term or explicitly introduced by hand to keep a count on the order of 
successive approximation involving powers of the perturbation-term Hj. In 
the latter case, it is set to unity, A — > 1, at the end of the computation to 
the desired order of approximation. 

In order that the above scheme works, at least the following obvious 
requirements must be met : 

(a) the system Hamiltonian permits a split in the above manner as in eqn.(l). 

(b) the perturbation-interaction term XHj remains sub-dominant to the 
solvable unperturbed part, H s . In view that the Hamiltonian is an operator, 
the above statement requires precise quantification. A commonly accepted 
criteria, which ensures the condition of sub-dominance of the perturbation 
term is the following 3 : 

< n s | H s | n s > | > | < n s \ XHj \ n s > | , (2) 

where | n s > is defined by the eigen value equation, 

H s \ n s > = E s n \n s > (3) 

(c) The method yields meaningful results if the actual physical spectrum 
can be computed as a power-series, i.e. 

E n = E s n + A EW + A 2 + , (4) 

and 

| n > = | n s > + A | n (1) > + A 2 | n (2) > + , (5) 

(d) It is further required that the above sequences represent either 
convergent or asymptotic series in order that the eqs.(4) and (5) make sense 
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as an approximation if truncation of the same after a finite number of terms 
would yield results representative of the actual physical situation. 

(e) An important corollary of the above requirement is, therefore, that the 
rate of convergence of the series ( or, the rate of decrease of the sub- 
asymptotic sequence, as the case may be ) are important considerations 
bearing upon the practical application of the method. 

The series-development in eqs.(4) and (5) can be generated by iteration start- 
ing from the Schrodinger equation for the states | n > and | n s >. Fol- 
lowing standard method 3 , one obtains the following general expression valid 
for arbitrary order : 

| n > = | n s > + Q n {H s -z)~ 1 (E n — XHj — z) \ n > , (6) 

After iteration, eqn.(2.6) yields the following result : 

I n > = Yl [Qn {H s -z)-\E n - \Hj - z) ] m | n s > . (7) 

m 

Similarly, for the energy shift one obtains : 
E n = E s n + ]T < n s \ XHj [Q n {H s -z)-\E n - XHj -z)\ m \n„> , (8) 

m 

where, z = arbitrary complex parameter; Q n = 1 — | n s >< n s |, 
is a projection operator with the property of projection to states orthogo- 
nal to | n s >, i.e. < n s \ Q n \ ip >= 0, for any arbitrary representative 
state \ ip > . The role of the arbitrary parameter 'z' is to generate 
by iteration, the two popular versions of perturbation expansion, namely the 
Rayleigh-Schrodinger(RS) 7 and the Brillouin-Wigner (BW) 10 versions for dif- 
ferent choice of z. Thus, for z = E n the RS-series is obtained while for 
z = E n , the BW-series follows. More explicitly, the separate expressions 
are given below for the two schemes : 

RS-Scheme 

\n > = E [ Qn (H s - E s n )-\E n - XHj 

m 

Similarly, for the energy-eigen value one obtains 
E n = E s n + £ < n s \ \Hj[Q n {H s -E s n Y\E n 



- E s n ) ] m I n s > , (9) 

- XHj - E s n ) } m I n s > , 

(10) 
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Similar expressions in the BW scheme are given below. 
BW-Scheme 

I n > = ]T [ Q n (H s - E n y\- \Hj) ] m \n s > (11) 

m 

E n = E s n + Y, < n *\ XR i [ Qn (H s - E n y\- XHj) ] m I n s > . (12) 

m 

It may be remarked here that the above (formal) expressions, which are 
valid to arbitrary order, become useful only when approximated to any given 
finite order. A comparision of the two schemes reveals that although the 
BW-expressions look simpler, yet it involves the unknown, exact energy of 
the system in the denominator, which thus requires further expansion up to 
the desired order in 'A'. This accounts for the popularity of the RS-scheme. 
On the other hand, however, there are systems for which, the BW-scheme is 
preferred as it converges more rapidly than the RS-scheme. 

Merits and limitations of perturbation theory 

Decidedly, the main merit of the method is the built-in provision for sys- 
tematic, order-by-order improvement over the zeroth-order (unperturbed) 
result. However, in its naive formulation as elaborated above ( which, we 
denote as naive perturbation theory (NPT) in the following) the applicability 
gets severly restricted due to several reasons. We mention some of these in 
the following. 

The main limitation arises when eqn.(2) fails to hold for the system. This 
may happen due to several reasons. Firstly, this may mean that the range 
of values of 'A' gets restricted (often, the condition: A < 1 ) provides a 
necessary but not sufficient condition. Secondly, the nature of the problem 
may be such that the inequality (eqn.(2)) fails ( even if A << 1 ). A third 
reason leading to the inapplicability of the NPT may arise when the desired 
condition, eqn.(2) does not hold for arbitrary value of the excitation label 'n'. 

Even when the basic requirement stated in eqn.(2) is satisfied, the appli- 
cability of the method may be restricted on the grounds of convergence of the 
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NPT. This may happen, for example, when the physical quantity sought to be 
calculated by perturbation series, is independently found to be non-analytic at 
A = 0. There is a large variety of physical phenomena, which fall under this 
catagory, e.g. tunneling, decay, phase-transitions, critical-phenomena, col- 
lective and cooperative phenomena such as super-conductivity, super-fluidity 
etc. It may be pointed-out in this context, that in most of the cases listed 
above, the perturbative ground state ( i.e., the ground state of the unper- 
turbed Hamiltonian, H s ) is shown to be unstable. In addition to the above 
limitations of the NPT on theoretical grounds, there may be practical difficul- 
ties arising, for example, from computing corrections beyond the first-order, 
which generally involve the (infinite) sum over the intermediate states. This 
problem has been addressed by several authors and it has been found in spe- 
cific examples that the difficulty can be surmounted by appealing to special 
techniques, such as the Dalgrano-Lewis method 11 . 

Some of the above aspects are exemplified when one considers the appli- 
cation of the NPT to the problem of the QAHO as discussed below. 

The QAHO problem and the NPT 

The Hamiltonian for the QAHO (in one space-dimension) is given by the 
following expression : 

H = l -p 2 + l -gx 2 + \x 2k ; k = 2, 3, 4, (13) 

In the above equation, p = and g, A are real, positive parameters. 

(The units are chosen such that m — 1 and H — 1). 

The NPT can be developed by recognizing that for A = 0, the system 
Hamiltonian reduces to that of the simple harmonic oscillator and is ex- 
actly solvable. Hence, it is natural to choose the anharmonic term as the 
perturbation-correction and express the Hamiltonian in the form given by 
equation (1), where XHj = \x 2k and H s = \p 2 + \gx 2 . 

In order to apply the NPT in the above formulation, it becomes apparent 
that the coupling strength be restricted to small values : A < 1 . However, 
it was discovered that the NPT-series fails to converge, even if the above 
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restriction is imposed. Specifically, it was demonstrated by the pioneering 
work by Bender and Wu 12 who studied the large-order behaviour of the NPT- 
series for the quartic-AHO and came to the conclusion that the co-efficient of 
X n grows as n! for any value of A even when A — > ! It was subsequently 
established 13 that the factorial-growth of the series worsened further for the 
cases of higher-anharmonicity. Even prior to the work under refs. [12,13], it 
was shown by Lam 14 that the development of the perturbation theory as a 
power-series in A has vanishing radius of convergence due to the occurence 
of an essential singularity at X = 0. 

Since then, a lot of investigation into the large-order behaviour of the NPT 
for the QAHO have been carried out and the results have found mention in 
texts 15 . From these investigations it has become clear that the divergence of 
the NPT could be ultimately traced to the failure of the essential condition 
expressed in eqn.(2) due to the eventual dominance of the perturbation-term 
Xx 2k over the unperturbed part for large amplitude of oscillation. 

In this context, it has now come to be recognized that the problem of 
divergence of the NPT may not be restricted to the case of the QAHO alone. 
For example ,in the case of quantum electrodynamics (QED), the possibility 
of failure of the NPT was conjectured by Dyson 16 much earlier. This has led 
to the conjecture that the problem of divergence of the NPT could be generic 
and the NPT-series may at most be an asymptotic series. A crude estimate 
of the optimal number, no of terms in the NPT-series beyond which, the 
subsequent terms start diverging could be obtained from the criteria : 

A no -R no « 0(1) , (14) 

where, it is assumed that the NPT-series (for the observable, say, energy E) 
is represented by 

E = J2 XnR m (15) 

n 

where, A is the coupling strength. Then, assuming factorial growth, R n ps n\ 
for large n, and using the Stirling-approximation, one obtains from eqn.(14), 
the desired estimate : 

n » 0(j). (16) 

It is clear, therefore, that the relevant coupling-strength must be sufficiently 
small in order that the NPT-series in eqn.(15) could make sense even as an 
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asymptotic, divergent series ! We elaborate on these aspects further in the 
following while discussing summability methods. ( It may be noted, in the 
context, that eqn.(16) may provide the rationale why the QED-perturbation 
series can still be trusted as an asymptotic series since the former starts di- 
verging only after ~ 137 terms since Xqed ~ 1/137 ). 

It may be relevant to note also the following aspects while discussing ap- 
plication of the NPT to the QAHO-problem : 

(i) For g = (see, eqn.(13)), the an-harmonic term : \x 2k can not be 
treated as a perturbation to the remaining part of the Hamiltonian, which 
is the kinetic-energy term (KET). This is because of the totally different na- 
ture of the resultant spectrum- the KET generates a continuous spectrum 
whereas the addition of the an-harmonic term changes this to a discrete one. 

(ii) Arguments based upon scaling and dimensional analysis 17 can be used for 
the QAHO-problem to show that the latter is essentially a single parametric- 
problem even though there are two couplings : g and A occuring in the 
Hamiltonian. Thus, for the quartic- AHO, the following scaling law holds 17 : 

E (g, A) = A§ E(g\-1, 1) , (17) 

i> (x ; g, A) = ij (x ^ ;g A~§, 1 ) , (18) 
where E = energy and ip = wave function of the system. 

(iii) The above scaling-property together with re-parameterization techniques 
can be effectively used 18 to develop perturbation theory for strong coupling 
regime, i.e., for A >> 1 since the effective single-coupling for the problem 
becomes g\~ 2 / 3 (see, eqns. (17-18). This is useful since the strong-coupling 
perturbation theory may be used as a complement to the usual small-coupling 
RS-series in studying smooth transition across A = 1. 

(iv) It may also be relevant to note that recently a novel-perturbation the- 
ory has been proposed by Turbiner 17 where the Schrodinger-equation is first 
converted to the Riccati-equation in the variable y{x) = ~ d ^f ^ > ( where 
ip(x) = wave-function) and then developing a perturbative expansion for 
y(x) through a recursive and iterative procedure. The resultant series for 
the energy and the wave function are characteristically different from those 
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obtained in the usual NPT and are claimed to be rapidly convergent. 



We close this sub-section noting that the difficulties encountered with the 
NPT as applied to the QAHO-problem have led to several alternative ap- 
proaches, as listed earlier which are discussed in subsequent sections. We 
next turn our attention to the method based upon variational-approximation. 

( 2 ) The Variational- Method of Approximation (VMA) 

The VMA is a powerful method in obtaining the approximate eigen-values 
and eigen-states of observables in QT. Like perturbation theory (PT), the ap- 
plication of the method has a long history- initial application of the method 
being ascribed to Lord Rayleigh 19 and to W.Ritz 20 . The VMA becomes espe- 
cially suitable when perturbation theory (PT) fails or becomes inapplicable 
to the considered problem. VMA is also used to test results based upon PT 
and for analysis of stability of the system considered. 

The essential ingredient of the method is based upon the variational theorem 
(Theorem-I) , which provides an upper bound for the ground state energy of 
the system : 

< ib I H lib > r , . , 

< W> s E w - E ° • (19) 

where, \ib > is any arbitrary, normalisable state; H is the Hamiltonian, 
and Eq is the ground-state energy of the system. ( In the above equation, 
we have used the standard Dirac-notation for expectation values and also 
denoted by E[ib], the energy- funtional). The proof of the theorem is based 
upon straight forward application of eigen-function expansion method in QT 
and can be found in any standard texts 3 . 



A second theorem (Theorem-II) found useful applications, states that the 
state | ib >, which renders the energy-functional E[ib } stationary i.e. 
SE [ib] = 0, satisfies the stationary Schrodinger equation : H\ ib > — E\ ib > 
Again the proof can be found in standard texts 3 and is therefore, omitted. 

For practical application to obtain a close approximation, one follws the 
Ritz-method 20 , which consists in specifying the trial -state | ib > as a 
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function of one or more free-parameters : {at} and then minimising the 
energy functional E[ip ] with respect to these parameters : 

dEtyydaii = ; d 2 E[iP}/da 2 > 0. (20) 



Thus, the Ritz-method enables one to obtain a least-upper bound (LUB), 
for a given choice of the trial-state, which is then the closest approximation 
to the ground-state energy for that choice. Some other aspects of the VMA 
are discussed below : 



(a) The equality sign in eqn.(19) holds if (and only if) the trial-state, | ip > 
coincides with the true ground state, \ ipo > of the system, (b) An error of 
I s * order, i.e., ~ 0(e) in the choice of the trial-state gives rise to an error 
of 2 nd order, i.e., ~ 0(e 2 ) in the computed value of the energy, which means 
that the energy is more accurately determined by the VMA than the wave 
function and consequently, the over-all accuracy of the method crucially de- 
pends on the ingenuity and insight in the choice of the trial-state, consistent 
with the physical boundary conditions for the system, (c) The generalization 
of the Rayleigh-Ritz method of the VMA to obtain estimates/approximation 
for the higher excited states can be achieved by several methods. Some of 
these methods are outlined below : 



(i) A tower of trial states, | if) a > ; a = 1, 2, 3 , approximating the 

higher excited states of the system may be constructed with the requirement 
that these are mutually orthogonal to each other, as well as, orthogonal to 
the variational-ground-state. The construction can be achieved by any suit- 
able method of orthogonalisation, such as the Schmidt-method 21 . Denote by 
E a , the energy- functional for the a-th state minimized with respect to its 
free-parameters and then arranged as an decreasing sequence, i.e., 



E 



a 



< 




H 









mm 



; a = 1, 2, 3, (21) 



The resulting sequence, E < E 1 < E 2 < E 3 , then represents the 

varitational approximation to the energy of the higher excited states in the 
order shown. 
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The other standard method to deal with the excited states within VMA 
is known as 22 the "method of linear variational approximation (LVA)". This 
consists of the choice of the trial-state as a linear superposition of a suitably 
chosen set of eigen-functions consistent with the boundary conditions and 
otherwise appropriate for the system, as follows : 

| i> > = 5Z C « I u n > 5 = 1 ' < u n \n m > = 5 nm . (22) 

n n 

The energy functional as given in eqn.(19) computed with this trial-state 
leads to the following set of equations after minimization with respect to the 
coefficients c n (see, eqn.(22)) : 

5>; (H nm - E5 nm ) = ; m = 1, 2, 3, (23) 

n 



Condition for existence of non-trivial solutions of the above set of equation, 
is given by the vanishing of the secular determinant : 

det ( H nm - ES nm ) = (24) 



For practical purpose, the determinant has to be truncated at some finite 
order l k\ The solution of eqn.(24) then provides, in general, /c-roots for the 
energy E, which then correspond to the approximate energies of the first 
l k — excited states'. It is obvious that the above method becomes efficient 
in practice when the truncation-error (due to the chosen finite value of k) 
becomes small and further that the spectrum stabilizes (with increasing value 
of k). Clearly, therefore, intelligent choice of the basis-states, | u n > becomes 
crucial in the above method. 

For certain problems, a variant of the above formalism is often called for 
when it becomes more convenient to work with a non-orthogonal basis set. 
In that case, eqn.(24) is replaced by the following equation : 

det ( H nm - ES nm ) = 0, (25) 

where, now S nrn = < u n \ u, m >, is the overlap function. 

The above basic formulation of VMA can be improved in several ways some 
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of which, are described in the following. 

(a) Improved choice of the trial- states 

This obvious generalization that has to be tried first for a given problem 
requires the choice of more realistic trial-states, which often involve large 
number of variation-parameters. These lead in general, to results that are 
more accurate. The main limitation of this approach is related to the fact 
that there is no systematic and standard method to follow. In the context of 
the QAHO, several trial-states have been tried with varying degree of suc- 
cess. A sample list is given under ref . [23] . 

(b) Obtaining lower bounds on the spectrum 

The idea in this approach is to obtain lower bounds to the energy eigen- 
values in addition to the upper bound provided by the standard VMA. Car- 
rying out the optimization of both the lower- and upper- bounds then leads 
more accurately to the actual eigen-value. Lower bounds are obtained by 
employing standard inequalities and other techniques. In the context of the 
QAHO some important work, in this direction, are listed in ref. [24] . 

(c) The Gaussian approximation 

This method employs the choice of a Gaussian as the trial state for a 
general problem. This approach relies upon the general applicability of the 
theory of small-oscillations to any arbitrary potential/interaction possessing 
a stable-minimum. The variation parameters are usually chosen to be the 
'width' and the location of the peak of Gaussian trial-state. For example, 
in one space-dimension, the trial wave-function for the ground state could 
be of the generic form, ip(x) = A exp[ ~^~ 2 J °^ ], where 1 xq and 'a' are 
the free-parameters and 'A J is a normalization constant. This method has a 
long-history mainly because it has been 'rediscovered' and 'refined' time-and- 
again. It has been applied to a veriety of problems in QT including quantum 
field theory. We refer to the work of Stevenson and collaborators 25 , which 
also provides guide to the earlier literature. 

(d) The method of minimum enrgy variance (MEV) 
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In stead of the Hamiltonian, H itself, upper bounds can be derived, in 
principle, for any arbitrary function f(H) by the generalization of the stan- 
dard Rayleigh-Ritz method. The generalization when the chosen function 
is considered as the variance of the Hamiltonian, has been termed 26 accord- 
ingly as "method of minimum energy variance (MEV)". The MEV has been 
shown 26 to result in improved estimate of the trial-state as well as, the energy 
eigen-value. In particular, it has been demonstrated that the estimates go 
beyond the simple Gaussian approximation 26 . 

In addition to the above, there are other methods, which combine the vari- 
ation method with other approximation methods, e.g., perturbation theory, 
super symmetry, the Hill-determinant method etc. in order to improve the 
results from the VMA alone. Some of these latter-schemes are separately 
discussed in the following subsections. 

( 3 ) The Hill- determinant and related methods 

In this method, the expectation value of the Hamiltonian of the prob- 
lem is computed in the trial state function, which is expanded in terms of 
of a suitably chosen set of basis functions (which need not be orthogonal) 
with unknown coefficients. The condition for the determination of these co- 
efficients translates into the condition of vanishing of the resulting secular 
determinant (Hill-determinant). This latter equation is then solved numer- 
ically for the approximants of the Hill-determinant obtained by truncation. 
More specifically, the trial-state proposed by the authors 27 , who pioneered 
this approach, is given as : 



ip(x) = exp (-x 2 /2)Y / c n x 2n . (26) 

n 

Substitution of this trial-state in the Schrddinger equation for the quartic- 
AHO then leads to the following three-term difference equation for the coef- 
ficients, c„'s : 



2 {n + 1) (2n + 1) c n _! + (E - 1 - An) c n - A c n _ 2 = 0, (27) 
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where, A is the an-harmonic coupling strength and E = energy of the system. 
( Alternatively, due to theorem-II mentioned earlier, the same result is ob- 
tained after minimisation of the energy functional, E [ip] with respect to the 
c* treating the latter as the variation parameters.) The condition of existence 
of solutions of the eqn.(27) then requires that the corresponding (infinite- 
dimensional) Hill-determinant D of the problem vanishes, D = 0.(The explicit 
form of D is given in ref. [27]). The authors of ref. [27] numerically solve the 
truncated approximant, D n = 0, by noting that these satisfy a three-term 
recursion relation : 

D n - (E-l-An) D n _i + 16An(n-l/2)(n-l)(n-3/2) J D n _ 3 = 0. (28) 

The solution, D n = 0, then leads, after ensuring convergence and stability, 
to n-roots for the energy E, which corresponds to the energy of the first n- 
levels of the even-parity states, when rearranged as an increasing sequence. 
( The energy levels of states with odd-parity can similarly be obtained by 
substituting, n — > n + 1/2). 

Obviously, therefore, stability, convergence, the computational time and 
effort needed, are the main criteria for the practical implementation of this 
method. The authors of ref. [27] obtained, by this method the eigen- values of 
the A :r 2fc -anharmonic oscillators separately for A; = 2, 3 and 4. The results 
compare well with the accurate numerical results obtained by other methods. 

However, several authors 28,29 have later pointed out the limitation of the 
approach in ref. [27]. One main problem concerns the normalizability of the 
resultant eigen-functions. Another problem is the failure of the method to re- 
produce the so called "terminated (polynomial) solutions" (later known as the 
"quasi-exact solutions" 30 ). Still another criticism 28 of the original method 27 
is the failure to reproduce correct spectral-properties when applied to a gen- 
eral an-harmonic oscillator having arbitrary couplings for the quadratic-, 
cubic- and quartic an-harmonic terms. Several modifications, involving im- 
proved choice of the variational trial-state, have been suggested 28,29 . For 
example, the following substitution of the exponential convergence factor 
(see, eqn.(26)), has been proposed : 

exp ( - x 2 /2 ) -> exp ( - a x 4 + (3 X 2 ), (29) 



16 



where, the constants a, (3 are either treated as variation parameters or 
determined otherwise. Even though the above 1 modified Hill-determinant 
method removes certain limitations of the original approach 27 , there still re- 
main certain controversy 31 regarding correct reproduction of the eigen-values 
and 'moments'. 

A related approach is the so-called "variational sturmain approximation 
" . The concept of the sturmain basis functions was introduced 32 ' 33 as a non- 
perturbative approximation scheme for solving the stationary Schrodinger 
equation in the context of several potentials considered in atomic- and molec- 
ular physics problems. This method involves, akin to the Hill-determinant 
method, the simultaneous solution of an infinite number of algebraic equation 
and, therefore, to the vanishing of the secular determinant as the necessary 
condition. For numerical solution it is necessary to solve the finite order ap- 
proximant to the secular determinant and thereafter ensure convergence and 
stability in respect of chosen size of the latter. Recently, it has been shown 34 
that the basic sturmain approximation can be variationally improved to yield 
better results. The authors in ref.[34] have obtained accurate results for the 
quartic-AHO and the exponential potential. For details, ref . [34] can be con- 
sulted, which also provides guide to the earlier literature. 

( 4 ) Combination of Variational and Perturbation techniques 

In order to retain the merits but to over come the limitations separately 
of perturbation theory and variation methods, several approaches have been 
proposed, where both the methods have been used in a single formulation. 
Various authors have worked using this combined approach, including the 
following: Halliday and Suranyi 35 , W.Caswell 36 , J. Killingbeck 37 , Hsue and 
Chern 38 , Feynman and Kleinert 39 and collaborators 40 , Patnaik 41 and Rath 42 . 
The review of each individual work as mentioned above falls beyond the 
scope of the present thesis. However, it may be relevant here to note certain 
common features of the techniques used by these authors and the consequent 
achievements attained. 

In most of the approaches listed above the basic Hamiltonian of the system 
is altered by the addition and subtraction of terms involving certain parame- 
ters. These parameters are apriori arbitrary but fixed later by the variational 
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minimisation of energy or by imposing other constraints, which simplify com- 
putation. Perturbation techniques can then be applied with redefinition of the 
unperturbed Hamiltonian and the perturbationcorrection. The resulting per- 
turbation series is often shown 35 ' 36 ' 41 ' 42 to be convergent. Thus, the problem 
of convergence of the naive perturbation theory as well as the absence of of 
a built-in mechanism of systematic improvement in the VMA, are overcome 
to a large extent in these hybrid method, which may be called variation- 
perturbation method (VPM). Apart from the above common underlying fea- 
ture, each approach differs in detail, which may be found in the individual 
reference cited above. 

(5) The operator Method of Approximation (OMA) 

The OMA has been pioneered by Feranchuk and Komarov 43 . The ap- 
plication of the OMA to the QAHO-problem and to several other problems 
has been described in ref.[44], which also provides a guide to earlier works 
of the authors. The basic idea of this method is to formulate the problem 
in the Fock-space of operators, instead of working with the co-ordinates and 
momenta. 

The relevant transformation is provided by the following relations : 

x = (a + a))/V2u ■ p = i{a ] - a)/yj{2/u), [ a, a f ] = 1 (30) 

where, a, a) are the annihilation- and creation-operators respectively and 
'u;' is a positive real number having the significance of the frequency of the 
associated simple-harmonic oscillator. The operators satisfy the standard 
commutation relation as shown in eqn.(30) in order to be compatible with 
the equal-time commutation relation, [x, p] = i. 

The development of the OMA then consists of the following steps : 

(i) The Hamiltonian of the system is first re-expressed in terms of the 
Fock- space operators ( by the help of eqn.(30)), 

(ii) The unperturbed Hamiltonian, H is chosen by clubbing all the diago- 
nal terms in H ( i.e., those involving polynomials of the number operator aJa), 
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(iii) The remaining off-diagonal terms are together chosen as the pertur- 
bation, V, 

(iv) A modified perturbation theory is then developed using the above 
choice of the unperturbed Hamiltonian and the perturbation term, 

(v) Finally, the arbitrary parameter 'u;' is fixed by variational minimiza- 
tion of the energy. The procedure is indicated as shown below : 

H(x,p) -> H(a,a!) = H (a^a) + V(a,(J) (31) 



To illustrate the method consider the QAHO problem with the Hamilto- 
nian given by 

H(x,p) = X -v 2 + \x 2 + A x 4 , A > 0. (32) 

Then the OMA-unperturbed Hamiltonian Hq and the perturbation term 
V, as defined in eqn.(31) are given by the following expression 44 : 

H = (1 + 2N)(uj 2 + l)/4w+ (3A/4cj 2 )(1 + 2iV + 2iV 2 ) , (33) 

V = -[(cu 2 -l)/4^]((a t ) 2 +a 2 ) + (A/4cu 2 )[2(a t ) 2 (2AT+3) + (2Ar+3)2(a 2 ) + (a t ) 4 +a 4 ] 

(34) 

In the above equations, N stands for the number operator, N = a)a. 
The energy eigen-value corresponding to Hq is easily obtained : 

= (n+l/2)[(^ 2 + l)/2^] + (3A/8cu 2 )(n 2 + n + l/2) (35) 

Minimization of the energy as given in eqn.(35) above then leads 44 to the 
following equation, for the frequency 'u;': 

to 3 - uj - 6A(1 + 2n + 2n 2 )/(l + 2n) = 0, (36) 

which determines this parameter as a function of the excitation level 'n' 
and the strength of the anharmonic coupling A. Substitution of the solu- 
tion of eqn.(36) in eqn.(35) then leads to the leading-order approximation of 
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the energy levels, which are found to be accurate within a few percent over 
the entire range of 'n' and A. Furthermore, this leading order result can be 
systematically improved order-by -order by developing the RS-perturbation 
theory with the perturbation term chosen to be V as given in eqn.(34). Be- 
sides, this perturbation series has been argued 44 to be convergent unlike the 
naive perturbation theory Similar results have been obtained 44 for the cases 
of higher an-harmonicity, double-well oscillators and a variety of other sys- 
tems. 

The main limitation of the model appears to be the implementation of 
the basic ansatz defining the split given in eqn.(31), particularly for non- 
polynomial interactions, interaction representing analytic functions express- 
ible in infinite series, problems in higher dimensions, quantum field the- 
ory etc. In such cases, additional assumptions/methods/skills have to be 
employed 44 . 

(6) The JWKB approximation method 

This method is quite old and popular due to its model-independent nature- 
it can be applied to any arbitrary, smooth potential. The formalism has been 
dealt at length in texts 45 . Therefore, it will not be repeated here. We note, 
however, the main features of the approximation scheme as applied to the 
discrete bound state problems in QT. The relevant central formula for the 
above purpose is the so-called JWKB quantization rule : 

/ k(x) dx = (n + -) 7T , n = 0, 1, 2, 3 (37) 

J a 2 

where, k(x) = 2m /h 2 (E — V(x)) ; a ,b are a set of adjacent "turning 
points" obtained by solving the equation, k(x) = and other terms have 
standard meaning. For those cases where the above integral can be evaluated, 
eqn.(37) can be expressed as : 

/(£, A, «?...) = (n + i) 7T, (38) 

where, the l.h.s. of eqn.(38) represents the value of the integral as a 
function of the energy (E) and coupling-strengths (A, <?....) occuring in the 
potential V(x). 
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The determination of energy E as a function of the other parameters 
(n, A, g...) therefore, requires the inversion of eqn.(38). This task is far from 
simple for most cases of interest since the inversion can not be achieved in 
closed, analytic form except in a few known cases ( which, also happen to be 
exactly solvable by standard methods ). However, eqn.(38) can be inverted 
by numerical-methods to the desired order of accuracy. The complexity of 
the procedure increases for the cases, for which the WKBJ-integral can not 
be obtained in closed form. Similarly, the cases involving multiple ( closely 
spaced ) turning points 46 and cases requiring higher order of approximation 
in the WKBJ-series, ( e.g. to deal with rapidly varying potentials) become 
increasingly difficult to implement. Nevertheless, several useful information 
/in sight about the energy levels can be obtained in limiting cases of the 
small-coupling regime, the strong- coupling and/or the large- n limits. Ap- 
plication to the case of the QAHO for studying the large- order behaviour 
of the naive perturbation theory was carried out in ref.[47]. The case of the 
DWO has been dealt in ref.[48]. The connection between multi-instanton 
solutions in the path-integral formalism, the JWKB method and large-order 
behaviour of perturbation theory has been discussed by several authors. We 
refer to Garg 46 who discusses the issue and lists earlier references. 

(7) Approximation methods based on quantum- canonical transformation 
(QCT) 

This is a powerful non-perturbative method, which has been successfully 
applied in QT, particularly in many-body systems exhibiting collective- and 
co-operative phenomena, e.g. superconductivity and super- fluidity 49-51 . The 
method had been primarily expounded by Bogoliubov 52 and hence more fa- 
miliarly known as the "Bogoliubov- Transformation" . These transformations 
connect the Hillbert-space of the interacting system to that of the interaction- 
free case, while preserving the canonical structure of the basic ( equal-time) 
commutation rules. The method is particularly useful in variational studies 
where the ansatz for the interacting vacuum state (IVS) generated through 
QCT can be tested to dynamically establish the stable ground state of the 
system. The basic formalism of the method can be illustrated, for a Bosonic- 
system with single-degree of freedom, through the following example. Con- 
sider such a system, which is described, in absence of interaction by the 
free-field annihilation- and creation operators, a,a^. The dynamic 'field' 
and the canonical momentum p are parameterized in terms of these oper- 
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ators as : = (a + at)/y/2 and p — % (at — a )/v / 2 and satisfy the 
standard equal-time commutation relation :[ , p ] = %. In presence of 
interaction, it is natural to assume that the interaction- free operators, a, at 
change to the interacting ones , b,tf . Similarly, let the vacuum states for 
the corresponding cases be denoted as : | > and \vac > respectively. By 
definition, these satisfy : a | >= 0, b \vac >= 0. The canonical- com- 
mutation rules satisfied by both the set of operators are identical by physical 
requirement, i.e., 

[ a , at ] = [ ft, &t ] = i (39) 

The general Bogoliubov-transformation relating the two sets of operators is 
given by 

b = a cosh(a) — at sinh(o;) , 

b^ = at cosh(a) — a sinh(o;) , (40) 

where, the parameter a is real and a priori arbitrary but the same can be 
determined dynamically through variational minimization of energy and/or 
by other physical requirements. The corresponding unitary transformation, 
which relates the two vacua is given by, 

| vac > = exp [1/2 tanh (a) (at at — a a) ] | > , (41) 

such that the parameter, a embodies the full effects of interaction : 
\vac > -> | > as a ^ 0. This vacuum state representing the situation 
in presence of interaction can be used as a trial state for a given Hamiltonian, 
such as given by eqn.(31) in order to study the stability and to compute the 
ground state energy. The energy for the excited-states can then be deter- 
mined by computing the expectation value: < n \H(a, at)\ n >, where the 
state | ra > is given by the standard formula : | n > = [(b^) n / \fnt] \v ac >. 

The formalism has been employed in refs. [53,54] for the QAHO and the 
DWO problems. In the context of A0 4 quantum-field theory, QCT has been 
employed in refs. [55,56,57]. QCT also provides important in-sight into the 
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structure of the interacting vacuum state 55 ' 56 ' 58 ' 59 . 

The limitations 60 of the method appear to be those of variational meth- 
ods as discussed earlier. 

(8) Approximation methods based on super symmetric quantum mechanics 

The method of super symmetric quantum mechanics ( SUSYQM ) has been 
dealt in texts and several review articles. As a representative text 6 , may be 
consulted, which provides guidance to earlier literature. As has been noted 
earlier, the property of shape-invariance in SUSYQM 6 has been successfully 
employed to considerably extend the class of exactly solvable potentials in 
QT. The method of SUSYQM can be profitably applied to improve and ex- 
tend the scope of known approximation methods, e.g. perturbation theory, 
variational method and JWKB-approximation scheme. In the context of 
non-relativistic QT, several exact results follow in SUSYQM valid for part- 
ner potentials 6 such as the property of iso-spectrality level degeneracy and 
positivity. These exact results provide a testing ground for approximation 
schemes applied to such potentials. 

The SUSY-improved perturbation theory starts from an initial guess of 
the ground state wave function, which can be based, for example, on a realis- 
tic variational ansatz. The super-potential corresponding to this trial-ground 
state wave function can then be constructed by computing the logarithmic 
derivative of the latter. The unperturbed Hamiltonian is then chosen to be 
the one obtained from the super potential and the perturbation correction is 
taken to be the difference of the original Hamiltonian and the unperturbed 
one. The development of the RS-perturbation theory then becomes straight- 
forword. For illustration of this method, ref.[6] can be consulted. 

The SUSY- based JWKB approximation provides the following modifica- 
tion of the quantization-rule : 

[ b \j2m [E { n ] - W 2 {x)} dx = nnh ; n = 0, 1, 2, 3,... (42) 

J a 

where, W is the super-potential and E$ is the n-th energy-level of the H±, 
which is one of the partner-Hamiltonians. Similar expression for the energy 
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levels of the other partner-Hamiltonian, H 2 is given by, 



/ \l2m [E ( n ] - W 2 (x)} dx = (n + l)nh ; n = 0, 1, 2, 3,... (43) 

It may be seen from the above equations that the exact result on the iso- 
spectrality / level-degeneracy of the partner potentials, is respected in the 
SUSY-JWKB quantization condition. Moreover, the formulae are valid for 
all values of 'n' rather than for large- n as in the case of the conventional 
formulation. The SUSY-JWKB formula is demonstrated 6 to yield the exact 
results for shape-invariant potentials and results with improved accuracy for 
other known cases. 

In the context of the QAHO/DWO, SUSYQM- methods have been inves- 
tigated by several authors. A sample list is provided in refs. [61-71]. 

(9) The methods based upon re-summation techniques 

The generic problem of divergence of the naive perturbation series (NPS) 
has led to the various re-summation techniques in order to extract some 
meaningful result out of the former. We will mention briefly the two popular 
re-summation techniques described in the literature, namely (i) the Borel-re- 
summation and (ii) the Pade- approximation in the following. 

(i) The Borel re-summation method 

The method of Borel re-summation 72 is a well-studied technique 73,74 , which 
permits, under suitable conditions, re-summation upto the optimal number 
of terms of an asymptotic power series exhibiting the growth at large-order 
'n' as n!. Consider such a series : 

F{g) = E fn 9 n , (44) 

n 

where, g is a real parameter ( 'coupling-strength' ) and f n ~ n\, for large- 
n. Then, the associated Borel re-summed series is defined as 

B{z) = J2fn* n /n\ , (45) 

n 
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This latter series is expected to have at least a finite radius of convergence 
by construction. The orginal series, (44) can be recovered, at least formally, 
from eqn.(44) as 

F(g) = (1/g) / exp (-z/g)B(z) dz , (46) 
J o 

by noting that / °° exp (—z/g)z n dz = n\g n+1 , provided that the integral 
on the r.h.s. of eqn.(46) can be uniquely evaluated. The integral represen- 
tation, eqn.(46) is known as the Borel-transform and this represents, within 
the radius of convergence of the series for B(z), the sum of the asymptotic 
perturbation series, eqn.(44). In order to be useful, the detailed knowledge 
of the analytic structure of B(z) in the complex- z (Borel-)plane is neces- 
sary. If the singularities of B(z) lie off the positive real axis then the integral, 
eqn.(45) can be evaluated by counter-integral after appropriate distortion of 
the contour around those singularities. 

However, if one ( or more) singularity in the Borel-plane lies on the real- 
axis, then the evaluation by contour integral becomes ambiguous; the integral 
representation, eqn.(46) no longer exists and, hence the case is not Borel re- 
summable. The other requirement is that the radius of convergence of B(z) 
must be larger than ( or equal to) the magnitude of the coupling-strength, 

in order that the re-summation is relevant. 

There is another situation, which invalidates the re-summation. This oc- 
curs when the terms in the original series all have equal phase. In that case 
the Borel-transform, B(z) itself grows too fast at infinity. Consequently, the 
Borel-integral, eqn.(46) can not be evaluated. There may be other situa- 
tions in addition to those already stated, when the singularity structure of 
B(z) prevents Borel re-summation. For the details, ref. [75] may be consulted. 

Finally, it may be stressed that the method fails when the coupling be- 
comes large, | g \ > 1 because then the original perturbation series becomes 
invalid due to the divergence of the individual terms. Thus, the re-summation 
techniques are necessarily restricted to the small coupling regime. 

Borel re-summability for the quartic-AHO has been proved in ref. [76]. 
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The relation of the singularity in the Borel-plane to the corresponding Eu- 
clidean classical action was first established by Lipatov 77 . An early study of 
Borel re-summation method was made by Jaffe 78 . Use of conformal mapping 
to enhance the rate of convergence of the Borel-transform has been discussed 
in ref. [79]. Bounds on the optical sum of divergent series, in the context 
of QT, have been obtained in ref. [80] using a novel variational method for 
the Borel-transform. Borel-summability of interacting quantum field theory 
has been reviewed, for example, in ref. [81]. Some modern applications of 
re-summation techniques have been discussed in ref. [82,83]. Singularities ( 
'Renormalons') in the coupling-plane due to the effects of renormalization of 
field theories and their significance have been discussed in ref. [84] . 

(ii) Pade- approximation 

Re-summation of the original perturbation series, eqn.(44) via Pade 
approximation can be expressed as the ratio of two polynomials, Pi(z) / Q m {z) 
of approximately same degree, as follows : 

[l/m\ F = Pi{z)/Q m {z) = (p +p 1 z+p 2 z 2 +...+p l z l )/(l+q 1 z+q 2 z 2 + ....+q m z m ) 

(47) 

The I + m + 1 coefficients, p , pi, p 2 , pr, q±, q2----<lm are chosen in 

such a way that the Taylor-series expansion of the r.h.s. of eqn.(46) agrees 
with the orginal perturbation expansion up to the desired order for small \z\, 
i.e., 

F(z) - [ Pi{z)/Q m (z) ] = O ( z l+m+1 ) ,z - (48) 



The condition, eqn.(48) requires the solution of l+m+l linear equations for 
the coefficients of the polynomials Pi(z) and Q m (z). If this system of equa- 
tions has a solution, the Pade-approximant can be expressed as a ratio of 

two determinants which depend upon the coefficients, / , /i, f 2 , fi+m+i 

of the original perturbation series, eqn.(44). The generalization to multi-point 
Pade-approximation is straightforward. If the original asymptotic expansion 
can be made about a set of given points, z = z a ; a = 0,1, 2,,.... 
then following the above procedure, the approximant to the original series 
can be developed for expansion about each point. Details may be found in 
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ref. [85,86] and references contained therein. 

Earlier works on Pade approximation for the AHO are contained in ref. [87, 
88]. It has been found 87 ' 88 that single-point Pade approximation works for 
the quartic- and sextic-AHO cases although the convergence is poor for the 
latter case. However, the same method fails 87-89 for the case of the octic- 
AHO. More recently, it has been shown in ref. [90], that clever use of scaling 
and transformation of the expansion variable results in accelerated conver- 
gence over a larger domain of z, by the use of two-point Pade approximants 
for the case of AHO's with quartic-, sextic- and octic an-harmonicity. It 
may be noted, however, that sophisticated numerical analysis using symbolic 
manipulation are required for obtaining results by the above method. 

(10) Approximation Schemes Based upon the Path-Integral Formulation of 
QT 

The Path-Integral (PI) formulation of Quantum Theory 91 " 93 has tradi- 
tionally served as an important basis to formulate and improve upon approx- 
imation methods- it may be recalled that the Feynman- diagram technique, 
which is perhaps the most popular method of computation in perturbation 
theory, originated from the PI- approach 94 . In the context of quantum me- 
chanics^ very comprehensive account of the PI- formulation is given in ref. [9] , 
which also provides extensive reference to earlier work. The approximation 
methods can be formulated either in the real-time formulation 91 ' 92 or as the 
imaginary time, Euclidean-formulation 9 . The latter is found to be ideally 
suited for extending the application to quantum-statistics. 

Various approximation schemes including the Semi-classical approximation 
methods 9 ' 95 , Stationary-phase approximation 9 ' 96 , loop-wise expansion 9 , large- 
order estimate of perturbation theory 9 ' 95 , variational estimation using 
inequalities 9 ' 96 , variational perturbation theory 9 and many more methods 9 
have been profitably based upon the Pi-formulation. Application to quan- 
tum field theory and statistical physics /critical phenomena are extensively 
treated in the text 95 . 

In the context of the QAHO-problem, the stationary-phase approxima- 
tion was used 97 in the Pi-formulation to infer the in-applicability of the naive 
perturbation theory due to the existence of an essential singularity at the ori- 
gin of the coupling strength plane. A variational method was proposed in 
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the Euclidean formalism of the PI in ref.[96] to obtain an approximation to 
the partition function of the QAHO and thereby, to obtain an estimate of 
the ground state energy. The method has been improved upon and extened 
in application in ref.[9] under the name of 'variational-perturbation theory 
(VPT)'. Wide applications and convergence tests of VPT are described in de- 
tail in ref.[9]. In ref. [98], variational lower-bound as well as, upper bound on 
energy of the QAHO have been obtained by the use of inequlities for the par- 
tition function. Semi-classical approximation of partition functions for one- 
dimensional potentials using the Pi-formalism has been described in ref. [99]. 
The same in higher dimensions is dealt in ref. [100]. The multi-instanton 
effects for potentials with degenerate minima is described in ref. [101]. In 
the context of scalar field theory, ref. [102] describes non-standard expansion 
techniques for the generating functional in the PI formulation. Several other 
applications in field theory and statistical physics can be found in ref. [95]. 

The major limitation of the Pi-formulation of approximation method appears 
to be the mathematical complexity involved in evaluating multi-dimensional 
integrals - even the simple, exactly solvable potentials such as the hydrogen 
atom problem, require rather complex techniques 9 for implementation. In 
many cases, it is therefore, preferable to use the standard operator based 
formalism for implementing the approximation methods. 

(11) The self- consistent schemes of approximation 

These schemes are popular in the context of many-body systems and other- 
wise known under the name of "mean-field approximation" and the "Hartree- 
approximation" including various generalizations of the latter, (for example, 
the "Hartree-Fock" method, " Hartree-Fock-Bogoliubov" method etc). The 
basic concept of the scheme consists in approximating the given potential 
by an exactly solvable one-body potential. If the original Hamiltonian is H 
and the approximating Hamiltonian is H , then perhaps the simplest way 
to ensure self-consistency in the approximation could be imposition of the 
following constraint : 

< n\H\n > = < n\H \n > , (49) 
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where, the states | n > are the eigen-states of Hq. In practice, the ap- 
proximating Hamiltonian may be chosen with a-priori unknown, adjustable 
parameters, which then get determined through the constraint, eqn.(49) and 
other requirements such as the variational principle and / or further simpli- 
fying conditions. The non-linear feedback that usually characterizes the self- 
consistency condition can then be easily ensured. Thus, the simple-looking 
equation, eqn.(49) possesses the potential to include interaction effects non- 
perturbatively and simultaneously preserve the non-linearity of the original 
Hamiltonian. The above method of implementation often yields results which 
are reasonably accurate even in the leading order of approximation, which 
can be further improved through standard recursive procedure as in the case 
of the Hartree-Scheme. 

In the context of A0 4 quantum field theory and the QAHO problem, such 
a scheme was considered in ref. [103] where, however, the expectation value 
(see,eqn.(49)) was restricted to the ground state only, in implementing the 
self-consistency requirement. The limitation has been overcome in proposing 
a "generalised Hartree-method" in ref. [58]. The 'mean-field' approach has 
also been tried in ref. [104] . 

The main limitations of the method appear to be the lack of uniqueness 
of the initial choice of the approximating Hamiltonian and the rate of con- 
vergence of the recursive procedure for subsequent improvement. 

(12) Other Methods of Approximation 

In addition to the methods listed above there are several other methods 
of approximation in QT. Since it is not possible to discuss all such methods 
within the limited scope of the present thesis, we list below only a small sam- 
ple which, we believe, are worthy of mention- the choice is entirely personal. 

(i) The 5 -expansion method 

This method 105 is similar in many respects to the approaches of the Gaussian 
approximation 25 , the method of Caswell 36 , Killingbeck 37 , Halliday and Suranyi 35 
and other methods discussed above under the category : ' Variation-perturbation 
methods 1 . The basic idea is to split the Hamiltonian into a 'free' part and 
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an 'interaction-part', both involving one ( or more) additional parameter, 
say Q, which is not present in the original Hamiltonian. In addition, an ex- 
pansion parameter, '<5' is multiplied to the interaction part for the purpose 
of development of perturbation theory. After calculation to the desired or- 
der, this parameter is set to unity, 5 — > 1 at the end. Hence, the method 
is intrinsically non-perturbative although calculation can be done order-by- 
order in the artificial parameter, S. The result of approximation to any finite 
order depends upon Q whereas the sum to all orders must be independent 
of the same. The sensitivity of the result to this parameter at any given 
order is minimized by application of the " principle of minimum sensitiv- 
ity (PMS)" 106 . The convergence of the 'perturbation theory' through the 
5-expansion method, has been rigorously established for the QAHO/QDWO 
problem 107 . Application to quantum field theories have been described in 
ref.[108]. However, the method, in its orginal formulation, fails to reproduce 
the wave-function at any order due to the variational nature of the approx- 
imation implied by the PMS. This defect has been remedied in ref.[109] by 
generalization of the scheme to allow co-ordinate dependence of the Q pa- 
rameter and using the PMS criteria also for the wave-function. 

(ii) The non-perturbative renormalization group(NPRG) method 

The NPRG formulation uses the Wilson-effective action (WEA) 110 as 
developed to study the critical phenomena, non-perturbative aspects in sta- 
tistical mechanics and quantum field theory. The method has been applied 
to quantum mechanics by Aoki and collaborators 111 . The basic equation in 
the formalism was derived by Wegner and Houghton (WH) 112 . However, the 
WH-equation, which represents the exact cut-off dependence of the WEA, 
can not be analytically solved exactly. Hence some approximation has to be 
made. In the so-called "local potential approximation (LPA)", quantum me- 
chanical systems, including the AHO and the DWO, have been numerically 
investigated 113 and shown to yield accurate results. Details can be found in 
ref.[113] and references to earlier related work contained therein. 

(iii) The method of multiple scale perturbation theory (MSPT) 

The above method has been propounded in ref.[114]. It is especially suitable 
to the cases where the ordinary perturbation theory diverges due to secular 
terms (terms that grow rapidly with the co-ordinate or the time variable) 
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in the chosen perturbation. The AHO /DWO problem is an obvious example 
where the perturbation Ax 4 soon predominates over the harmonic term, gx 2 
no matter how small A is. MSPT recognizes different characteristic proper- 
ties of the system, such as the frequency of oscillation, different asymptotic 
behaviour of the wave function etc, at different time/length-scales when sec- 
ular terms are present. In ref. [114], the MSPT techniques have been applied 
to solve the Heisenberg-equations of motion for the QAHO problem, obtain 
the energy spectrum and the eigen-states. It also provides references to re- 
lated earlier work. 

(iv) The coupled cluster expansion method (CCEM) 

The CCEM is a known technique in many-body theory. In that area, the 
CCEM has been widely applied 115 and extensively reviewed 116 . The basic 
idea is due to Hubbard 117 who proposed that the true interacting vacuum 
state of the system can be generated by an operator, exp(S) acting on the 
free, non-interacting vacuum state. In that respect, the above relation can be 
generated through appropriate canonical transformation, as discussed above. 
The method has been extended to quantum mechanics and quantum field 
theory in ref. [118]. In ref. [119], the technique has been applied to a general 
AHO with quadratic, cubic and quartic self-interaction such that various 
cases of QHO/DWO are covered. The starting point is a parameterization 
of the interacting vacuum, given by the relation : 

\vac' > = exp ( so) + to) 2 )|0 > , (50) 

where, s and t are real parameters and other terms have standard meaning 
as explained earlier. The creation/annihilation operators b, for the inter- 
acting vacuum state \vad > can be constructed by standard method 119 . 
Using the above trial vacuum state and the tower of higher excited states 
generated there from, the usual variational calculation can be done 119 for the 
AHO and it yields fairly accurate result as discussed earlier. 

For the situation of degenerate ground state as happens in case of the DWO, 
the trial vacuum-state can be constructed as a linear superposition of the 
degenerate vacua and variational calculation can be made to calculate the 
energy split as well as the average energy of a given level. The CCEM method 
consists in systematically improving the above initial trial state through a 
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series of linked cluster operators parameterized by the relation : 

oo 

\vac > = exp ( ^2 Sk)\vac > , (51) 
k = i 

where, Sk = Ukip^) k . Computation of the energy eigen values of quantum 
systems can be based upon the above trial state and higher excited states 
based upon it after suitable approximation to restrict the infinite sum in the 
argument of the expoential in eqn.(51). For details, ref.[119] can be consulted. 

This completes our survey of the various approximation methods in QT. 

In summary, it is revealed from the above brief yet hopefully representa- 
tive survey that considerable progress has been achieved in the development 
of approximation schemes to tackle various problems in quantum theory. At 
the same time it is also brought to focus by the above survey that the 'ideal ' 
method of approximation, based upon the earlier stated, desirable set of re- 
quirements of general applicability, simplicity, systematic improvability and 
efficiency of computation etc, remains largely elusive till date. In view of the 
above situation, there still remains a lot of scope to incorporate, if possible, 
the above set of desirable criteria into a single-scheme. 

The present dissertation is devoted to an attempt at achieving the above 
ambitious objective, namely, to develop a self-consistent, non-perturbative, 
simple yet efficient scheme of approximation, which is systematically improv- 
able and generally applicable, in principle to an arbitrary interacting system 
in QT including the quantum field theory. The scope of the present dis- 
sertation is, however, limited in application to the self-interacting Bosonic- 
systems in quantum theory in which, we include the quartic an-harmonic 
oscillator, quartic double-well potential, sextic an-harmonic and double-well 
oscillator, the octic-an-harmonic oscillator and the A0 4 quantum field theory 
without spontaneous symmetry breaking. In the next Chapter, we describe 
the formulation 120 . 



3. The New General Approximation Scheme 
(NGAS) 
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3.1 Formulation of NGAS for self- interacting Systems 

Consider a generic Hamilitonian H\(4>,p) describing a self-interacting 
quantum system involving the field '0' and conjugate momentum 'p', given by 

H x (<f>,p) = H s (^p) + \H I ( ( j ) ), (52) 

where, H s is the unperturbed Hamiltonian and XHj is the self- interaction 
with A as the coupling strength. We use the language of field theory identi- 
fying quantum mechanics as field theory in (0 + 1) dimensions. 

Many physically important systems are described by the above Hamilto- 
nian including the anharmonic oscillators ( AHO ); the double- well oscillators 
( DWO ); the A0 4 quantum field theory in the symmetric-phase as well as, 
in the spontaneously symmetry broken (SSB) phase ; pure Yang-Mills fields 
with the quartic ' gluori- self interaction etc. 

As has been discussed earlier, perturbation theory treating H s (<f),p) as 
the unperturbed Hamiltonian and the entire interaction \Hi(<f)) as pertur- 
bation, fails to converge even for infinitesimal A . Other non-perturbative 
methods often suffer from limitations as reviewed in Chapter-2. With a 
view to overcome these limitations and difficulties the formulation of NGAS 
consists of the following steps: 

3.1.1 Choice of the approximating potential : V(4>) 

The aim is to replace the orginal interaction \Hi(4>) by a "suitably cho- 
sen" approximating potential \V(<f)) such that it preserves the symmetries 
of the original interaction but leads to an " exactly solvable" Hamiltonian, 
i.e. the "Effective Hamiltonian (EH)" generated by V((f>) and defined by 

H (<j>,p) = H s (<f>,p) + \V(<t>), (53) 

is exactly solvable, i.e. 

H \n> = E^ ] \n>, < m\n > = 5 mn , (54) 

where, the spectrum \n > and the eigen- values E^\ are known. ( We 
consider, for simplicity, that the spectrum is discrete and non-degenerate ). 
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We refer to this requirement, eqn.(54), as the "condition of exact solvability" 
(CES). The next requirement on V((f>) is the condition of equal quantum 
average, as explained below. 

3.1.2 The Principle of Equal Quantum Average 

The effective Hamilitonian is constrained to yield the same quantum av- 
erage (QA) as the original, i.e. 

< n\H((j>)\n> = < n\H ((f))\n >, (55) 

which implies, 

<n\H I ((f))\n> = <n\V((j))\n> (56) 

We refer to this requirement, eqs.( 55-56) as the "condition of equal quan- 
tum average ( CEQA)". The next step is to optimize the approximation as 
described below: 

3.1.3 Variational Optimisation of V((j>) 

Let V{4>) involve a set of free-parameters: { }. Then the requirement 
of optimisation consists in the variational minimisation of H with respect 
to the free-parameters (cej) which characterize V((f>) : 

A < H„ >= (57) 

where, the notation is 

< A > = < n\A\n > (58) 
We refer to this condition, eqn.(57) as the "condition of optimality (CO)". 

3.1.4 The Self-Consistency and other features of the Method 

The steps which are outlined above are the essential ingredients of the 
proposed approximation scheme in the leading order (LO). The following sev- 
eral observations regarding the approximation scheme are in order: 
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(i) It is to be noted that in a restricted form, i.e. when the quantum average 
in eqn.(57) is restricted to the ground state only, the CEQA, as expressed in 
eqn.(55), corresponds to the Hartree-approximation/ mean field approxima- 
tion in quantum field theory 103 . In view of this, the NGAS can be regarded 
as a "generalised" Hartree-approximation method 58 . 

(ii) The self- consistency of the procedure is implicit in eqs.( 53-55): the 
states | n > which are obtained as the solution of Hq (see, eqn.(54) ), are 
used as input (in eqn.(56)) to determine V{4>) which, in turn, defines H 
(eqn.53), thus making the "feed-back loop" complete. This can be schemat- 
ically represented as: 



(iii) The leading order (LO) approximation in NGAS consists of finding the 
spectrum \n > and the energy eigen-values E®. This is easily achieved 
because of the CES (eqn.(54)). It may be emphasized, however, that even the 
LO results capture the dominant contribution of the full (nonlinear) interac- 
tion through the requirement, eqn.(55), eventhough one always deals with an 
exactly solvable Hamiltonian, H . We consider this as a key feature of the 
approximation method. 

3.1.5 The Vanishing Quantum Average of the Effective 
Interaction and its Significance 

The other important result which follows trivially from eqn.(55) is that 
the modified interaction XH' defined by the relation: 



In otherwords, only the non-diagonal matrix elements of H' are non- van- 
ishing in the basis which renders H Q diagonal. 

This result is of considerable significance and it naturally suggests a 
scheme of improved perturbation theory (IPT) in which, H is chosen as 




XH' = X(H T -V), 
has vanishing QA for arbitrary 'A' and L n\ i.e. 



(59) 



< n\XH'\n >= 0. 



(60) 
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the unperturbed Hamiltonian and \H' is considered to be the perturba- 
tion. The IPT thus developed, can be shown to be convergent, owing mainly 
to eqn.(60) and thus can be used to systematically improve the LO result 
order-by-order (Chapter 10). The various steps outlined in equations (53- 
60) define the scheme in the LO. 

In the following Chapters, we implement and illustrate the general ap- 
proach described above to specific quantum systems with self-interaction. 



4. Application of NGAS to the case of the 
Quart ic-Anharmonic Oscillator: formulation for 
the leading order(LO) 

4.1 Application of the NGAS to the case of the Quar- 
tic Anharmonic Oscillator in Quantum Theory 



In the following subsections we illustrate the NGAS applied to the case 
of the quartic anharmonic oscillator(QAHO) in quantum theory. However, it 
may be worth while to first describe the importance of the QAHO problem 
due to its applications in diverse areas of quantum physics and its signifi- 
cance from a theoretical stand-point as well. 

4.1.1 The Quartic Anharmonic Oscillator(QAHO) in Quantum The- 
ory: its importance and applications 

The quantisation of the simple harmonic oscillator is perhaps among the 
most celebrated examples 8,9 due to its immense significance and application. 
In particular, small-oscillations about the equlibrium- position are described 
by (normal-mode) simple-harmonic oscillators. Thus, quantum theory of 
the simple harmonic oscillator describes small oscillations of complex dy- 
namical systems about the equlibrium. This property leads, for example, 
to quantisation of free-field systems, which are shown equivalent to an as- 
sembly of independent harmonic oscillators. Similarly, most lattice-dynamic 
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studies in the area of condensed matter physics employ the harmonic ap- 
proximation as the leading description of the underlying physics. However, 
realistic approximation of actual physical systems need the incorporation 
of the anharmonic-interaction when deviation from small-oscillations or in- 
clusion of self-interaction among oscillators become necessary on physical 
ground. This is the primary reason which endows the study of anharmonic 
oscillators with considerable importance and significance and explains its 
application in diverse areas of physics 8,9 . Anharmonic oscillators are, there- 
fore, among the most widely studied models in quantum theory in which 
the quartic-anharmonic oscillator ( QAHO ) is the simplest system exhibit- 
ing self-interaction. This system has been extensively investigated leading 
to a vast amount of literature 8,9 . Its importance arises due to physical ap- 
plications in areas including field theory, condensed matter physics 121 , sta- 
tistical mechanics 122 , non-linear systems 123 , classical and quantum chaos 124 , 
inflationary cosmology 125 , lattice dynamics, plasma oscillations etc, to cite 
only a few cases. Besides, the QAHO has also served as a theoretical lab- 
oratory to study convergence of perturbation theory 126 , development of 
non-perturbative approximation methods 127 , renormalisation 25 ' 113 , vacuum 
structure 128 and stability analysis 129 etc. 

The Hamiltonian of the system is given by 

H = \ V 2 + I 9 2 + A0 4 , (61) 

where A , g are real and positive. Note that the canonical momentum 
conjugate to the 'field' <fi(t) is given by p(t) = d<j)(t)/dt . In the notation of 
eqn.(52), the free field Hamiltonian corresponds to 

H s {4>) = \p 2 + \g 2 (62) 

and the interaction is 

A# 7 (0) = A0 4 . (63) 

In order to develop the NGAS for the QAHO, we follow the steps outlined 
in the Chapter-3, as described below. 

4.2 Formulation of the NGAS for the QAHO 
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4.2.1 Choice of V(</>) 



The first step is to find a suitable approximating potential V((f>). The 
following ansatz is naturally suggested on the grounds of simplicity and exact 
solvability: 

V ((f)) = A (j) 2 - B <f) + C. (64) 

where A,B,C are parameters to be determined self consistently. It would 
appear from eqn.(64) that the global symmetry of the orginal Hamiltonian, 
eqn.(61), under <fi — > —(f), is not respected by the ansatz, eqn.(64). However, 
this is illusory since the co-efficient 'B' in eqn.(64) is dynamically determined 
to be propertional to < (f> > ( see eqn. (87) below), which ensures the preser- 
vation of the above global symmetry. 

4.2.2 Solution of the eigen value equation of the Resultant 
Effective Hamiltonian(EH) 



To obtain the exact analytic solution of the spectrum of V(<f>) we pro- 
ceed as follows. 

We define the "Effective Hamiltonian (EH)" by the following equation 

#o(<M = \v 2 + ^0 2 + AV(0) (65) 

The EH defined in eqn. (65) can now be recast in the form of a diagonalizable 
structure as outlined below. Substitution of eqn. (64) in eqn. (65) leads to the 
following equation 

#o(<M = \v 2 + I 9 2 + HA 2 - B + C) (66) 
The simplification of the eqn. (66) leads into the following Hamiltonian: 

H = + I w 2 (0 - a) 2 + Zio, (67) 

where, 

uj 2 = g + 2\A, (68) 
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a = XB/uj 2 , (69) 

and 

ho = \C - W. (70) 

It may be atonce recognized that the EH given by eqn.(67) corresponds to 
a "shifted", effective, harmonic oscillator where both the field, as well as, 
the energy are respectively shifted by 'a' and '/io'. Note further that the 
c-number parameters 'u;' and '<x' are restricted by physical requirement, to 
satisfy uj > 0; er = real (since, = Diagonalisation of iJ can be 
achieved by the standard method of invoking the creation-and annihilation 
operators, defined by 

4>{t) = a + (b + 6+)/>/2ZJ, (71) 

p(0=<v^/2(6 t -6), (72) 

along with the standard equal-time (canonical) commutation relation (ETCR) 
given by 

[ M f ] = 1- (73) 

We first introduce the number operator: = tfb and its eigen-states by 
the equation 

N b \n > = n\n> . (74) 

Because of the hermiticity and positivity of the operator N b ; the eigen values 
are positive-definite and the eigen states are orthonormal: 

< m \n > = 5 mn , (75) 
^|nxn| = I. (76) 

n 

Further, as proven in standard text books 130 , by the use of the ETCR (eqn. 
73), the spectrum of eigen values is restricted to: 

n = 0,1,2, (77) 

Then, using the defining parametrisation in eqns. (71-72) one immedi- 
ately obtains the Hamiltonian in diagonal form: 
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H = cu (N b + 1/2) + ho. (78) 
The energy spectrum is then trivially obtained and given by 

= iv^ + h , (79) 
where £ = (n + 1/2) with 'n' taking values according to eqn. (77). 

At this stage, the following remarks/ observations and results may be 
noted: 

(i) By using the standard properties of the creation-/annihilation opertaors, 
the QA of polynomials of field '0' and momentum 'p' can be easily evaluated. 
In particular, we note the following results for subsequent use (in the follow- 
ing, the vacuum-state is denoted by \vac > and is defined by the property, 
b\vac > = 0): 

< vac\(f)\vac > = < n\<f>\n > = < >= a , (80) 

< 02 > = a 2 + ^ /uj) . <p 2 > = ^ (gl) 

< 3 > = cr 3 + 3 <t(£/uj), (82) 

< 4 > = a 4 + 6 a 2 (e/cu) + 3 (l + 4e 2 )/8cu 2 , (83) 

<06> = a 6 + 15a\Z/u) + 45a 2 (l + 4Z 2 )/8u 2 

+ (5/8)(e/c 3 )(5 + 4e 2 ). (84) 

The QA of the original Hamiltonian H defined in eqn. (61) can now be 
evaluated with the help of eqs. ( 80-84) and given by 

<H> = \<p 2 > +^#<0 2 >+A<0 4 >. (85) 

Substituting the QA values from eqs. (80, 81, 83) and noting that the 
eqn. (55) holds, eqn. (85) is transformed to: 

<n\H\n> = <n\H \n>=u£/2 + (g + l2\a 2 )(£/2Lu) 

+ (3A/8cu 2 )(l + 4£ 2 ) + ga 2 /2 + Xa\ (86) 
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(ii) Eqn.(80) shows that '<r' corresponds to the vacuum expectation value 
(VEV) of '0'. In view of this result and eqn.(69), the coefficient 'B' of the 
linear term in V(<f>) (see eqn.(64)) can be re-expressed as : 

B = (uj 2 /X) < > . (87) 

Eqn.(87) demonstrates that the global symmetry of the original Hamil- 
tonian under the transformation, — > — 0, is preserved by the poten- 
tial V((f>) , which is not otherwise transparent in eqn.(64). 

(iii) If one denotes by a and a' , the corresponding creation- and annihi- 
lation operators of the 'free' theory (i.e. defined by H s (<f>,p)), then '0' and 
'p' can also be expressed in terms of these operators analogous to eqs. (71-72) 
as: 

0(t) = a+ (a + a^/v 7 ^, (88) 

p(t) = i^/2 (a f -a), (89) 

where ojq = y/g . It is important to note here that both the sets of 
creation- and annihilation operators satisfy identical (equal- time) commu- 
tation relation: 

[a, a f ] = 1 = [6, &t] (90) 

It may also be mentioned that the VEV of (< > = a) remains 
invariant in the two descriptions. Eqn.(90) implies that the two sets of oper- 
ators must be related by a quantum canonical transformation ( " Bogoliubov 
transformation") 52 . This result has crucial implications for the vacuum struc- 
ture and stability of the approximate theory. This is discussed in detail in 
Chapter 9. 

Returning to the implementation of the NGAS, the next task is to deter- 
mine the free parameters, 'a' and 'cj' (or equivalently, A, B and C occuring 
in eqn.(64)). This is achieved as follows: 

4.2.3 Determination of the free parameters 

The conditions expressed under CEQA and CO given in eqs( 56-57) are 
sufficient to fully determine the free parameters involved in the approxima- 
tion. These requirements translate in this case, to the following equations: 
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< 4 > = A < <p 2 > -B < > +C ; (91) 

d < H > /did = and d < H > /da = where, < H > 
is given by eqn.(86). Carrying out the explicit variational minimisation of 
< Hq > with respect to 'a/ i.e., d < H > /duj — one obtains the following 
equation : 

cu 3 - u(l2\a 2 + g) - 6A/(£) = 0, (92) 

Again carrying out the explicit variational minimisation of < H > with 
respect to 'cr' i.e., d < H > /da = , the following equation is obtained: 

a(4Aa 2 + g + 12A^/cj) = 0, (93) 

where, in eqn. (92), /(^) = £ + l/4£. Eqs. (92-93) are to be solved simultane- 
ously to determine 'cu' and '<t' as functions of A, g and £. In the following, we 
refer to these eqs. (92-93) as the "gap equation (GE)" and " the equation for 
the ground state (EGS )" respectively. The constants A, B and C appearing 
in eqn. (64) can now be determined as follows: 

To determine the parameter 'A' we use the eqs. (68) and (92). Rearranging 
the "gap-equation" given in eqn(92) we obtain the relation 

^ = ((? + 12Aa 2 ) + ^^ (94) 

Eqn. (94), considered together with eqn. (68) lead to a determination of 
the constant 'A', given by 

A= 6 a 2 + 3/(0 /u (95) 

To determine the parameter 'B' we use the eqs. (69) and (93). Rearranging 
the eqn. (69) leads to the following relation: 

<r(-l + ^) = (96) 
Comparing eqs. (93) and (96) we obtain the parameter 'B' as 

B = (1 + g) (ato 2 /X) + 4uj 2 a 3 + 12w<rf (97) 



42 



To obtain the parameter £ C we use equation(91): 

C =< 4 > -A <(f) 2 > +B <</)>, (98) 

On substitution of the values < <fi >, < <f) 2 >, and < 4 > 
from eqs.(80, 81, 83), along with the values of 'A' and 'B' determined above, 
the parameter 'C is determined. Using the gap-equation (92) further the 
parameter C is given by 

-6A/(0 + 12A£ - ^4^] - <x 4 (5 + 8c 2 ) (99) 
On substitution of eqn.(99) in eqn.(70) we obtain the relation 



3A , . 2 \ 3A£/(£) o r , a; 



-6A/«) + 12\( - - A<7 4 (5 + fc 2 ) (100) 

Thus, the parameters defining the approximating potential V((f>) are com- 
pletely determined. 

Next, using eqn.(86) and noting that < H > = < H > = one 
obtains the energy spectrum: 



=^/2 + (( ? +12A ( T 2 )(e/2c) + (3A/8c 2 )(l+4a+^ 2 /2 + Aa 4 , (101) 

where 'uj^ and 'cr' are solutions of eqs. (92-93). The different solutions thus 
obtained are discussed in the following: 

4.2.4 The Leading-Order (LO) Results - Determination of the 
Spectrum of H Q 

Solution of the gap-equation (GE) eqn.(92) and the equation for ground 
state (EGS) eqn.(93) constitute the key ingredients in the calculation of the 
energy spectrum. It is convenient to first obtain the solution of the EGS, 
eqn.(93). For the case of the QAHO ( g, A > ) , eqn.(93) has two solutions: 

(i) (j = (102) 
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and 



(ii) g + 4Acr 2 + 12\£/u) = 0. (103) 

We now analyse the eqn.(103) which gives 

a 2 = -^{g+l2K/uj) < (104) 

Since g, A > 0, eqn.(104) does not lead to acceptable solution, which 
requires a 2 > . Therefore, the only physically accepted solution for 
QAHO is the solution (i) given in eqn.(102), i.e., a = 0. ( This is also 
intuitively obvious since the single-well shape of the "classical" potential can 
not get altered to double- well shape by quantum - fluctuations). Substitution 
of eqn.(102) in eqn.(92), then leads to the following simplified GE for the 
QAHO: 

^ 3 -^-6A/(O=0. (105) 

It may be emphasized at this point that this GE ( with g=l ) has been de- 
rived by several authors 131 , but starting from widely different considerations. 
On substitution of a = in eqn.(101) we obtain the equation 

*r -^+i(f)+|> + 4a. doe) 

Using the GE, eqn.(105) and noting that 

/(0 = i± ^f 1 , (107) 
the last term in eqn.(106) can be simplified to be 



8^1 +4a = (|)(^) (108) 

Substitution of eqn.(108) in eqn.(106) then leads to the following simple 
expression for the energy spectrum of the QAHO in the LO: 

^ = |(3o; + ^) (109) 

where 'to' is obtained as a solution of the GE for the QAHO given by 
eqn.(105). It is to be noted that the "gap-equation" given by eqn.(105) is in 
the form of a cubic equation of the type 
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x 3 - 3Px -2Q = 0; P, Q > (110) 
The real solution of this eqn.(HO) is given by, 

x = Q 1/3 [(i + ^(i- J)) 1/3 + (i - ^(i-^)) 1/3 ] (in) 

Then comparing the coefficients of eqs.(HO) and (105) the real solution 
is obtained explicitly as 

= (3A/(0) 1/3 [( 1 + Va-P)) 173 + ( 1 - \f0^ r p)) 1/ % (H2) 

where, p~ 4 = 243A 2 / ' 2 (£,)/g 3 ■ It may be noted that, for the case 
when g = 1 , the solution for ' uo ' as given above has the correct limiting 
behaviour, u — ► 1 for A — > and further that it exhibits the non 
- analytic dependence on the coupling A at the orgin characteristic of the 
non-perturbative nature of the NGAS. 

As has been noted by several authors 131 , the formula, eqn.(109) is ac- 
curate to within a few percent of the 'exact' result for the full allowed range 
of g, A > (both in the 'weak coupling-' and the 'strong coupling' regimes) 
and for all values of the excitaion level n > 0. In particular, the accuracy 
in the strong coupling regime can be judged by the following result for the 
computation of the ground state energy. The 'exact' asymptotic result ( for 
g = 1) is given by 44 , for A — > oo, E \ exact = 0.668A 1 / 3 which is to be 
compared to the LO-result in NGAS: Eq \ ngas = 0.681A 1//3 in the same limit. 
In Table-1, the LO results of the present approximation scheme for g = 1 are 
presented for sample values of 'A' and 'n' along with 'exact' results ( obtained 
by numerical methods (Hsue and Chern, ref.[23])) and results from earlier 
computaion 44 ( obtained by different analytical methods) for comparision. 
As can be seen from the comparison, the LO results, lie within 0.2 — 2 % of 
the results from the ' exact numerical calculations 23 . In the same Table-1 , 
we also dispaly the improvement of results obtained by inclusion of the first 
non-trivial correction in the improved perturbation theory (IPT), which is 
discussed in Chapter-10. 

We next apply the method to the case of the quartic -double well oscillator 
in the following Chapter. 
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Table- 1: Leading order (LO) results and perturbation correction in the 2 nd 
order of IPT for the energy levels of the QAHO computed for g = 1 and 
sample values of 'A' and 'n' shown along with analogous results of ref.[44], 
compared with the 'exact '(numerical) results of (Hsue and Chern ref.[23j). 
The relative percentage errors in the two schemes are also shown. 



A 


n 




Exact 


E^ 


Error(%) 


£(2) 
n 


Error(%) 








ref.[23] 






ref.[44] 


ref.[44] 


0.1 





0.5603 


0.5591 


0.5591 


0.007 


0.5591 


0.007 




1 


1.7734 


1.7695 


1.7694 


0.005 


1.7694 


0.005 




2 


3.1382 


3.1386 


3.1391 


0.016 


2.9006 


7.580 




4 


6.2052 


6.2203 


6.2239 


0.058 


5.4795 


11.96 




10 


17.266 


17.352 


17.374 


0.127 


14.539 


16.32 




40 


94.843 


90.562 


95.766 


5.75 


76.152 


15.91 






0.8125 


0.8038 


0.8032 


0.070 


0.8032 


0.07 


1 


2.7599 


2.7379 


2.7367 


0.043 


2.7367 


0.043 


2 


5.1724 


5.1792 


5.1824 


0.061 


4.4440 


14.19 


4 


10.900 


10.964 


10.982 


0.17 


8.8890 


18.93 


10 


32.663 


32.933 


33.013 


0.243 


25.833 


21.56 


40 


192.79 


194.60 


195.15 


0.282 


149.87 


22.99 






1.5313 


1.5050 


1.5030 





.131 


1.5030 


0.131 


1 


5.3821 


5.3216 


5.3177 





.070 


5.3177 


0.070 


2 


10.324 


10.347 


10.356 





.090 


8.6131 


16.76 


4 


22.248 


22.409 


22.457 





.210 


17.651 


21.23 


10 


68.171 


68.804 


68.996 





.270 


52.943 


23.05 


40 


409.89 


413.94 


415.18 


0. 


.300 


316.13 


23.62 



100.0 


3.1924 


3.1314 


3.1266 


0.150 


3.1266 


0.150 


1 


11.325 


11.187 


11.178 


0.080 


11.178 


0.080 


2 


21.853 


21.907 


21.927 


0.090 


18.095 


17.40 


4 


47.349 


47.707 


47.817 


0.230 


37.314 


21.80 


10 


145.84 


147.23 


147.65 


0.285 


112.79 


23.40 


40 


880.55 


889.32 


892.03 
46 


0.300 


677.91 


23.77 



5. Application to the case of the Quartic 
Double Well Oscillator (QDWO): 
Leading Order(LO) Results 



5.1 Importance of the QDWO and its Applications 

The QDWO is also an extensively studied system 132 because of its theo- 
retical importance and practical application. The Hamiltonian of the system 
is given by: 

H = l -p 2 - \g<? + \<\>\ g,X >0 (113) 

The crucial -ve sign of the (f> 2 - term generates a quite different physical 
situation than the case of the QAHO, even in the classical limit. The 'clas- 
sical' potential, V c = — \g<t> 2 + A0 4 exhibits the familiar double- well shape 
with symmetric minima. For g — 1 , these minima are located at positions 
± — ^ and each with depth ^ ■ ^ s ^ becomes smaller and smaller 
the depth of the two- wells become deeper and deeper. The actual low lying 
energy eigen-states of the problem become radically different from the trial 
wave function offered by the harmonic basis. This fact severely handicaps the 
convergence of the resultant perturbation theory. In principle, the natural 
solution would be the simultaneous use of two harmonic basis centers around 
the positions of the minima at ± . However, the implementation of 
that idea, although possible, implies the use of nonorthogonal states which 
is rather cumbersome. 



The other difficulty is that the theory is not defined for A — > , be- 
cause the ground state does not exist in that limit due to the non-existence 
of a lower limit to V c . In that sense, the SHO is not the free-field limit of 
the QDWO. Therefore, the naive perturbation theory (NPT) is not appli- 
cable as such, to this case. The perturbation expansion of the eigen values 
E n (g,X) in powers of A is divergent 12 ' 13 for all A > 0. This fact may be 
understood qualitatively by noting that the addition of the term Ax 4 turns 
a completely continuous eigen- value spectrum of p 2 — g x 2 into a com- 



47 



pletely discrete spectrum bounded from below. A nonperturbative treatmant 
is therefore, necessary. The WKB method is one such method which is well 
suited especially for lower energy eigen values. However, as shown below, 
the NGAS can be sucessfully applied to the case of the QDWO in completely 
analogous procedure as in the case of the QAHO, for a considerable larger 
range of values of n,g and A. 

5.2 Application of NGAS to the QDWO 

5.2.1 The Gap Equation for the QDWO 

To develop the NGAS for QDWO we use the identical ansatz for V(4>) 
as given in eqn.(64). The parameters A, B, C can also be determined self- 
consistently as before. The " Effective Hamiltonian(EH)" for the system is 
analogously defined as: 

#oOM = \v 2 - \g<P 2 + W{<P)- A, g > (114) 

Using the ansatz for V(<f)) given in eqn.(64) in eqn.(114), the EH for 
QDWO is then recast into the diagonal structure analogous to the case of 
QAHO and is given by 

Ho = \ v 2 + (0 - a) 2 + ho, (115) 

where, the parameters are defined as before: 

u 2 = g + 2\A, (116) 
a = XB/uj 2 , (117) 

ho = AC - -uj 2 a 2 . (118) 

In this case also, the EH defined in eqn.(115) is to be interpreted as a 
"shifted" harmonic oscillator. Also in this case, the frequency uj > ; and 
the vacuum configuration-parameter a = real, are the physical require- 
ment. 
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Using eqs. (71,72,73) and proceeding as in the case of the QAHO, one can 
express H into the desired diagonal form given by the eqn.(78) from which 
one obtains the energy spectrum 

Eg* = ut + ho, (119) 
where f = (ra + 1/2) ; n = 0,1,2, 

In order to implement the NGAS, the next task is the determination of the 
free parameters 'cr' and 'a;' (or, equivalently, A, B and C occuring in eqn.(64) 
for the QDWO. 

For this purpose, taking the quantum average of original Hamiltonian 
defined in eqn.(113) and using eqs. (80) (81) (83) leads the following equation 

<n\H\n> = <n\H \n>=uj£/2- ^ g (a 2 + 

+A(a 4 + 6 a 2 (£/u) + 3(1+4 £ 2 )/8u 2 ), (120) 

Carrying out the variational minimisation of EH given in eqn.(120) with 
respect to free parameter 'u;' i.e., d < H > /duj = leads to the following 
"gap-equation (GE)" for the QDWO: 

u 3 -u(12\ ( r 2 -g)-6\f(O = 0, (121) 

Similar procedure with respect to the free parameter 'a' i.e., d < H > 
/da = , leads to the following equation for the ground state(EGS) con- 
figuration: 

a(AXa 2 -g+ 12A£/cu) = 0. (122) 

It is here to be noted that these equations differ from the analogous equa- 
tions, (92-93) for the QAHO by the substitution, g — > —g, as expected. 

It is again convenient to solve the EGS first. It should be noted however 
that, in contrast to the case of the QAHO, there are now two realizable 
quantum phases of the system corresponding to the solution of eqn.(122) for 
the ground states. These are characterised by the solutions: 

4Aa 2 = £ - 12A(£M (123) 
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and 



a = 0. (124) 

respectively. The two solutions lead to different "quantum phases" for 
the QDWO, which is discussed below. 

5.2.2 Different Quantum- Phases of the QDWO and the 
Critical Coupling 

Following standard terminology used in the literature, the solution 
given by eqn.(123) leads to the " Spontaneously Symmetry Broken (SSB)" 
phase whereas, the other solution, eqn.(124), corresponds to the " Symmetry- 
Restored ( SR )" phase. It is shown below that the dynamic realization of 
the two "phases" is controlled by the coupling 'A' such that the SSB phase 
is energetically favoured when A < A c , whereas the SR phase is preferred 
for A > A c where, A c is a ' critical' coupling. To demonstrate this we consider 
the GE in the respective phases: 

5.2.3 Solution of the Gap Equation and Determination of the 
Spectrum in the Different Phases 

(i) The SSB-phase of the QDWO 
The GE for SSB-phase is obtained by substitution of eqn.(123) in eqn.(121) 
and is given by 

ul - 2gu a + 6Ap(0 = 0, (125) 

where, = 5£ — l/(4£) and we have denoted by u a , the frequency 
in the SSB phase. To get physicallly accepted solution of eqn.(125), it is 
convenient to adopt the 'trial' solution given by: 

u a = p [e 10 ' 3 + e~ ie ' 3 } = 2p cos(0/3). (126) 

where 'p' and are to be determined. Substitution of eqn.(126) in 
eqn.(125) leads to a determination of 'p' and cos 9 : 

p = y/(2g/3), (127) 

and 
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cos 9 = - (A/A c ), (128) 

with 

X c =(2g/3f/ 2 /3p(0 = X c (g,0 (129) 
From eqn.(128), one gets either: 

9 = n + cos-\\/\ c ), (130) 

or 

9 = 7 ^ + sin- 1 (X/\ c ). (131) 

However, on stability grounds it is found that eqn.(131) leads to the 
acceptable solution as the corresponding ground state energy lies lower than 
that obtained using the other solution, eqn.(130). Thus the physical solution 
of eqn.(125) is given by 

^ = 2^|)co S [^ + i S m- 1 (A)] (132) 

An estimate of 'A c ' for the ground state and for g — 1 is A c (g — 1, £ = 1/2) 
= 0.0362886. Clearly, the solution given in eqn.(131) is valid only when, 
X < A c . 



The energy-levels in the LO of the SSB-Phase are computed in analogous 
manner to the QAHO as follows. One notes that 



e {o)\ Q dwo = <n \H^ DWO \n> = <n\H$ DWO \n> 



+ o- 2 {-\g+^) + \o-\ (133) 

Z CO 

where H$ DWO is given by eqn.(114) and the QA values, as given in 
eqs. (80-83) have been substituted. From the defining eqn. of the SSB - 
Phase i.e., eqn. (123) one can substitute for a 2 and obtain: 

^l? S T°4^ + f + ||(l-20a-^ (134) 
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The third-term in eqn.(134) can be further simplified by the use of the GE 
in eqn.(125) and given by 

»L(1 = <!>(,«.-£). (135) 

Substitution of eqn.(135) in eqn.(134) then leads to the final expression for 
the energy levels given by 

E^r= (|)( 3 ^+^)-(^)- ( i36 ) 

(ii) The SR- Phase of the QDWO 

For A > A c , the SR-Phase is dynamically favoured. The stability anal- 
ysis of the phase-structure confirms this. The GE, in this case, is obtained 
by substituing eqn.(124) in eqn.(121): 

u 3 s +gu s -6\f(O = 0, (137) 

where, we have denoted by u s , the frequency in the SR-Phase. Note that 
the above equation simply follows from the GE of the QAHO, eqn. (105), by 
the substitution: g — > —g, as expected, due to the underlying single well 
shape. The energy levels in this phase can be calculated rather easily. On 
substitution of eqn. (124) in eqn. (133) one obtains the equation: 

F (0)\QDWO _ n\H QDWO \n > I 

\SR = < n \ n \ n > l<7=0 

= 5"-f-5# + S (1 + 4{1) (138) 

Rearranging the GE for SR-Phase given in eqn. (137) leads to the following 
equation for the last term in eqn. (138): 

!^(l + 4a = (|)(^ + -)- (139) 

QLOg 4 UJ S 

Using the result obtained in eqn. (139) in the simplified eqn. (138), the energy 
levels for SR-Phase are given by the following simple expression: 

EW Q sS WO = (kto.-±). (140) 

4 Us, 
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which, again follows from the corresponding formula for the QAHO, 
eqn.(109), by the substitution, g — > —g . In eqn.(140), oo s is the solu- 
tion of eqn.(137) which can be obtained in analogous manner and is given 
by 

u s = (3A/(0) 1/3 [( V(l+P) + l) 1/3 " (y/(l+P)~ 1) 1/3 ] (141) 

where, p' 1 = 243A 2 / 2 (0/# 3 . 
In Table-2, we present the energy-levels of the QDWO in the LO, over a 
wide range of 'A' and 'n' for g — 1. The results are compared with an earlier 
computation 36 , which employs a modified perturbation theory and includes 
correction up to twenty orders of perturbation. In the same Table-2, we 
also dispaly the improvement of results obtained by inclusion of the first, 
non-trivial correction in the improved perturbation theory (IPT), which is 
discussed in Chapter 10. From the comparison with earlier calculation 36 , 
it is seen that the LO results are already quite accurate. 

In the next Chapter we consider the case of the sextic- anharmonic os- 
cillator and the sextic- double- well oscillator in the NGAS. 
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Table- 2: The computed energy levels of the quartic - DWO in the lowest 
order of NGAS for sample values of '\ ' and 'n' compared with the results 
of ref.[36] which includes perturbation correction up to twenty orders in a 
"modified" perturbation theory. Also shown are the results obtained after 
inclusion of the perturbation correction at the next order in IPT. 
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n 




E^ 


ref.[36] 


n 1 

U . J. 
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U.1UUU 
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n 77fn 

U.I 1 uo 
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1 ^fi^fi 

X .ouou 


1 (^47 


1 .uouu 




4 


o.OoUD 


Q 70Q0 


o.OoUZ 




10 


12.192 


12.517 


12.400 


1.0 





0.5989 


0.5752 


0.5800 




1 


2.1250 


2.0800 


2.1800 




2 


4.2324 


4.2600 


4.2500 




4 


9.4680 


9.5950 


9.5600 




10 


30.530 


30.650 


30.420 


10.0 





1.4098 


1.3752 


1.3800 




1 


5.0650 


4.9910 


5.0900 




2 


9.8660 


9.9050 


9.8900 




4 


21.561 


21.791 


21.700 




10 


66.950 


67.820 


67.620 


100.0 





3.1340 


3.0650 


3.0700 




1 


11.175 


11.024 


11.002 




2 


21.638 


21.715 


21.700 




4 


47.023 


47.505 


47.200 




10 


145.27 


147.10 


146.70 
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6. Application of NGAS to the Sextic-Anharmonic 
and Double Well Oscillator 



6.1 Importance of the Sextic-Anharmonic Oscillator 
and its Applications 

The sextic AHO-system is an example of higher anharmonicity, which 
is also widely investigated 12,13,133 . This system is interesting and impor- 
tant in its own right as it finds application in diverse areas of physics. As 
a theoretical laboratory, this system together with its double-well counter- 
part provide perhaps the simplest examples, for which supersymmetric quan- 
tum mechanics (SUSYQM) has definite predictions 61-71 for energy-levels for 
specific values of the Hamiltonian-parameters. Hence, various models and 
approximation methods can be tested against the exact analytic results of 
SUSYQM. Besides, owing to the higher anharmonicity, the divergence of the 
naive-perturbation theory at higher-orders becomes even more severe 133 , 
E n ~ T[n(m — 1)] for \<f) 2m type of AHO. This result endows the sextic- 
AHO with added importance for testing convergent-approximation methods. 
It may also be noted that the "Wick-ordering" method of Caswell 36 is not di- 
rectly applicable to this case since the method generates 4 counter-terms 
which are not present in the original (non-ordered) Hamiltonian. For the 
above stated reasons of practical applicability and theoretical importance, 
the sextic-AHO provides a unique testing ground for the NGAS, which is 
described in the following. 



6.2 Application of NGAS to the Sextic-AHO 
6.2.1 The Gap Equation for the Sextic-AHO 

The Hamiltonian for the system is given by 

H = \v 2 + I 9 4> 2 + A0 6 , (142) 
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where, A ,g > 0. As in other cases, 

H s (<t>) = \p 2 + \g 2 , (143) 

corresponds to the free-field Hamiltonian. The interaction-term is, there- 
fore, 

\Hr((/>) = A0 6 . (144) 

To apply the NGAS, we follow identical steps as in previous cases. An 
identical ansatz is assumed for the potential V((f>) , which is given in 
eqn.(64). Hence the Effective Hamiltonian (EH) for this case is analogously 
defined as: 

#o = \v 2 + I 9 4? + A V{<f>) (145) 

Substituting the ansatz for V(<f>) in eqn.(64) in eqn.(145), the EH can 
be recast into the diagonizable structure, as in earlier cases: 

H = \p 2 + \uj 2 (0-a) 2 + h , (146) 
where, analogously 

lu 2 = g + 2\A, (147) 



a = XB/cu 2 , (148) 

ho = \C - \uo 2 a 2 . (149) 

The interpretation of Hq in eqn.(146) is also identical- it is the Hamilto- 
nian for a "shifted" effective harmonic oscillator. Following steps identical to 
those in the case of quartic-oscillators, the diagonal form for H is obtained 
as: 

Ho = to (N b + 1/2) + ho , (150) 

where again N b = tfb with b, W defined as before (see eqs. (71-72)). The 
eigen value equation as given in eqn.(54) then leads to the energy spectrum 
as given below: 

= ut + ho, (151) 
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where £ = (n + 1/2); n = 0,1,2, and h is given by eqn.(149). 

For the determination of the parameters 'u;' and 'a', we use, as before, the 
variational minimisation conditions: d < Hq > /du = and d < Hq > 
I da = . For this purpose the quantum average of eqn.(142) evaluated as: 



<H> =<H > = \<p 2 > + \ g < 2 > +A < 6 > (152) 

On substitution of the results given in eqs. (80,81,84) in eqn.(152), the 
following equation results: 



< n\H\n > = < n\H \n>= u£/2 + ^ g (a 2 + 

Carrying out variational minimisation of eqn.(153) with respect to the 
frequency parameter 'u;' we obtain the relation 



u 4 -cu 2 (# + 30Aa 4 )-45A(crV2O(l + 4£ 2 )-(15A/4)(5 + 4£ 2 ) = 0, (154) 

Further, carrying out the variational minimisation with respect to £ <j' 
leads the following equation: 



10£a 2 15(4£ 2 +1) N1 
a[g + 6A(cr + + ' ) ] = 0. (155) 



Eqn.(155) which defines 'a' provides the corresponding ground state con- 
figuration(EGS) of this system. 

As in earlier cases, it is convenient to first analyse the EGS which leads 
to the following relations: 

(j = (156) 

and 

g + 6A(cr 4 + 10fr 2 /u + 15(4£ 2 + 1)/8uj 2 ) = 0. (157) 
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On inspection of eqn.(157), it is easily seen that no physical solution for 
a 2 (i.e. a 2 > ) exists for g > 0, A > 0. Hence, eqn.(157) is not 
physically realisable for the sextic AHO. Therefore, the physical ground state 
of the sextic AHO is uniquely determined by the l a = 0' solution as given by 
eqn.(156). 

6.2.2 Solution of the Gap Equation and Determination of the 
Energy Spectrum 

The GE given in eqn.(154) simplifies to the following form on substitution 
of the correct EGS expression, i.e. a = : 

c/-£cj 2 - (15A/4)(5 + 4£ 2 ) = (158) 

The parameters A, B, C defining V(<f)) can be determined in analo- 
gous manner as given below. Rearranging the gap eqn.(154) and using the 
eqn.(147) 'A' is dtermined as 

A= 15 tr 4 + 45<r 2 (l + 4£ 2 )/4£cc> + (15/8cj 2 )(5 + 4£ 2 ); (159) 

Similarly, following the same procedure, 'B' is calculated using eqn(148) 
and the EGS given in eqn.(155) and given by 

B = a[(l + g)(co 2 /X) + 6cjV + 60a 2 u^ + (45/4) (1 + 4£ 2 ) ]; (160) 
To determine 'C we use the equation as given below: 

C= < 6 > -A < (f) 2 > +B < >, (161) 

together with expressions for < 6 >, < <fi 2 > and < > as given in 
eqs. (80-84). On substitution of the physical solution, eqn.(156) for the EGS 
in eqn.(153) and using the GE as given in eqn.(158) one obtains the following 
simple expression for the energy-levels of the sextic AHO in the LO: 

E^\ sexU ^AHo = |(2w + -); g> 0. (162) 

where, ofcourse, l uJ , is the solution of the GE given in eqn.(158). In 
Table-3, sample results for the energy levels of the sextic AHO computed in 
LO are presented for g — 1, over a wide range of 'A' and 'n' and compared to 
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the results 133 ( shown within parenthesis) with percentage deviation from the 
latter (shown within square bracket). It can be seen from this tabulation that 
the LO- results obtained in the NGAS are quite accurate, compared to those 
obtained by sophisticated numerical calculations 133 . Further improvement in 
accuracy is achieved by application of the IPT as discussed in Chapter-10. 
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Table- 3: Sample results in the lowest order(LO) of (NGAS) for the sextic 
- AHO over a wide range of '\ ' and 'n' compared for g = 1 with the 
results of ref. [133] (shown in the parentheses). The relative percentage error 
is shown in square brackets. 
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A = 0.2 
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10.0 





1.193 


1.676 


2.323 




(1.174) 
[1.611] 


(1.610) 
[4.079] 


(2.206) 
[5.313] 


1 


3.966 


5.931 


8.420 




(3.901) 
[1.681] 


(5.749) 
[3.165] 


(8.115) 
[3.762] 


2 


7.420 


11.61 


16.74 




(7.382) 
[0.523] 


(11.54) 
[0.612] 


(16.64) 
[0.6179] 


4 


16.15 


26.48 


38.73 




(16.30) 
[0.9170] 


(26.83) 
[1.302] 


(39.29) 
[1.426] 


6 


26.88 


45.08 


66.36 




(27.29) 
[1.50] 


(45.94) 
[1.870] 


(67.70) 
[1.98] 


10 


53.24 


91.17 


135.0 




(54.31) 
[1.967] 


(93.26) 
[2.245] 


(138.2) 
[2.323] 


14 


85.01 


147.0 


218.3 




(86.78) 
[2.047] 


(150.4) 
[2.230] 


(223.4) 
[2.279] 


17 


111.9 


194.4 


289.0 




(114.0) 
[1.868] 


(198.3) 
[1.974] 


(294.19 
[2.001] 
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14.52 


20.39 


30.37 


(13.95) 


(19.56) 


(29.12) 


[4.148] 


[4.244] 


[4.298] 


29.16 


41.03 


61.18 


(28.98) 


(40.78) 


(60.81) 


[0.6157] 


[0.6145] 


[0.6138] 


68.01 


95.90 


143.2 


(69.05) 


(97.38) 


(145.4) 


[1.499] 


[1.517] 


[1.527] 


117.0 


165.1 


246.5 


(119.4) 


(168.5) 


(251.7) 


[2.043] 


[2.058] 


[2.067] 


238.7 


337.1 


503.8 


(244.5) 


(345.3) 


(516.1) 


[2.367] 


[2.377] 


[2.383] 


386.6 


546.2 


816.3 


(395.7) 


(559.1) 


(835.6) 


[2.306] 


[2.313] 


[2.316] 


512.1 


723.7 


1082.0 


(522.7) 


(738.6) 


(1104.0) 


[2.016] 


[2.020] 


[2.022] 



6.3 The Sextic-Double Well Oscillator (Sextic-DWO) 



6.3.1 Importance of Sextic-Double Well Oscillator 

The sextic double well oscillator has been the subject of much investiga- 
tion and discussion for a pretty long time 134 , both from the analytical and 
the numerical points of view. This is because it has important applications 
in quantum field theory and molecular physics 135 and in other branches also. 
By using Hill-determinant method, Biswas et al 27 have calculated the eigen 
values of the oscillators of the type \x 2n . Banerjee and Bhattacharjee have 
obtained energy eigenvalues for the potential of the form x 2 + Xx A using 
scaled Hill determinant method 136 . However, it has been pointed out that 13 
this method has limited domain of applicability for the sextic oscillator po- 
tential of the type V(x) = fix 2 — \x 4 + r/x 6 ; rj > . There have 
been some improvements 138 of the procedure for removal 139 ' 140 of dificulties 
in the Hill-determinant approach 141 . The potential cited above is particu- 
larly of great interest in scalar field theory 142 and also in the calculations 
of the vibrational spectra of molecules 143 . It has been discussed by authors 
of ref. [144] that the potential admits exact analytic solutions for the ground 
state under certain conditions. It has also been shown that 145 the poten- 
tial is "quasi exactly solvable (QES)" i.e., exact solutions can be obtained if 
the coupling constants satisfy some constraints. Moreover, supersymmetric 
quantum mechanics (SUSYQM) can be applied 71 to compute the eigen val- 
ues of the double-well potential x 6 — 3a; 2 and the spectrum can be related 
to that of the "partner" potential: x 6 + 3x 2 , which is the anharmonic 
counter-part . We apply our scheme of approximation in the next section, 
to calculate the energy eigen values for the sextic-DWO as a method of con- 
firmation of our scheme, as well as, to test the various exact results following 
from QES and SUSYQM. 

6.3.2 Application of NGAS to the Sextic-DWO 

The Hamiltonian for such a system is given by the following expression: 

H = \ V 2 - I 9 2 + A0 6 , (163) 
where, A > and g > . To apply NGAS, we follow identical steps 
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as in previous case of the QDWO. An identical ansatz is suggested for the 
potential V((f>) which is given in eqn.(64). Hence the ( Effective Hamiltonian 
(EH)) for this case can be written as: 

H = \ P 2 - \g 2 + A V{4>) (164) 

On Substitution the value for V(<p) given in eqn.(64) in eqn.(164) the EH 
is now expressed into the diagonizable structure: 

Ho = \ v 2 + \ (0 - <t) 2 + ho, (165) 

where, 

oo 2 = -g + 2\A, (166) 

a = XB/oo 2 , (167) 

h = XC - \uj 2 <j 2 . (168) 

The parameters A,B,C can be determined through identical procedure 
as in earlier cases. In particular, we get back relations given by eqs. (159-161) 
with the substitution: g < — > — g. Again the eqn.(165) can be regarded 
as the hamiltonian for the "shifted" effective harmonic oscillator for sextic- 
DWO, which is identical with the equations developed earlier for the case 
of QAHO eqn.(67), QDWO (eqn.( 115) and sextic-AHO (eqn.(146)). This, 
therefore, confirms the genraral applicability of the method formulated in 
Chapter 3. As usual ' oo ' is identified as the frequency of the "shifted" 
harmonic oscillator which has the restriction of satisfying the physical re- 
quirement oo > . The parameter ' a ' accounts for the field-shift of the 
"shifted" harmonic sextic-DWO which must be real, since the field is hermi- 
tian, = 0^. Parameters ' u ' and ' a ' are determined by following the 
analogous procedure ( see below ). Finally the EH defined in eqn.(165) is 
transformed by standard method into diagonalised form as: 

Ho = oo (N b + 1/2) + ho , (169) 

where the creation- and annihilation operators and ETCR given in eqs. (71- 
73) have been used. The energy spectrum for sextic-DWO is then trivially 
obtained as: 
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4°) = u£ + ho, (170) 
where £ = (n + 1/2) ; n = 0, 1,2, and h is given by eqn.(168). 

(i) The Gap-Equation and the different Quantum- Phases of the 
Sextic- QDWO 



For the determination of the parameter ' uo ' and ' a ' we use the stan- 
dard variational-minimisation conditions d < Hq > /duo = and 
d < H > I do = . On taking the quantum average of eqn.(163) and 
taking the advantage of eqn.(55) we obtain < H > as given below: 

< n\H\n > = < n\H \n >= uC,/2 - - g (a 2 + (£ / 'u)) 

+A(a 6 + 15 a' + 45 a\l + 4 £ 2 )/8c 2 + 5$(4 £ 3 +5) ), (171) 

From the variational minimisation of eqn.(171) with respect to ' to ', we 
obtain the relation 



-cj 2 (-^ + 30A ( t 4 )-45A( ( jV2O(1 + 4£ 2 )-(15A/4)(5 + 4£ 2 ) =0, (172) 

Eqn.(172) corresponds to "gap-equation" GE for the case of sextic-DWO. 
Similarly, from the variational minimisation of eqn.(171) with respect to ' o ' 
leads to the following equation: 

a \-g + 6A(ct 4 + I0£a 2 /u + 15(4£ 2 + l)/8u 2 ) ] = 0. (173) 

Eqn.(173) has the significance of defining the ground state (EGS). This equa- 
tion leads, as in the case of the QDWO, to the SSB-phase and the 
SR-phase, corresponding to the 'a 2 ^ 0' and the '<r = 0' solutions respec- 
tively. It is found on computation, however, that the SR- phase is energeti- 
cally favoured for all values of X since the energy levels for this phase, always 
lie below the corresponding ones for the SSB- phase. This result has impor- 
tant consequence in the context of supersymmetry, as discussed subsequently. 
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(ii) Solution of the Gap-Equation and determination of the en- 
ergy spectrum of the SR-phase 



The GE, in the case of sextic-QDWO for the SR-phase, is simply obtained 
by substituting l cr = 0' in the "gap-equation" eqn.(172) which corresponds 
to the GE for the AHO (eqn.(158)) by the replacement, g — > — g : 

^ 4 + ^-(^)(5 + 4£ 2 )=0, (174) 

where again, we have distinguished the frequency of the DWO, by the 
subscript 'a'. To determine the energy levels one substitutes ' a = 0' in the 
eqn.(171) for < H > one obtains: 

E^Uuc-dwo = <H >\ a=0 = ^ - i(^) + p|(5 + 4e 2 ). (175) 

Z Z UJ a SCO* 

Next, using the GE, eqn.(174), the last term in eqn.(175) can be simplified 

as: 

8c< 6 uj a 

Substitution of the above equation in eqn(175) then leads to energy spectrum 
given below: 

E^\ sextlc - DWO = |(2cu a -—),£> 0. (177) 

In Table-4, sample results for the energy levels of the sextic DWO com- 
puted in the LO are presented and are compared to the results of ref.[71]. It 
can be seen from this tabulation that the LO- results obtained in NGAS are 
quite accurate, compared to those obtained by sophisticated numerical calcu- 
lations of ref.[71]. On comparision it is seen that there is good agreement for 
all values of 'n', expect for the ground state. This is further discussed in the 
next Chapter. Further improvement in accuracy is achieved by application 
of IPT as discussed in Chapter-10. 

In the Chapter-7, we consider the implication of supper symmetric quan- 
tum mechanics (SUSYQM) for the case of the sextic oscillators and compar- 
ision with the results of NGAS. 
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Table- 4: Sample results for the energy levels of the sextic- AHO and DWO 
in the LO of NGAS (for (3 = 1) displaying the approximate validity of ISPP 
relation(see eqn.(191)). Also shown for comparision, are the corresponding 
results of ref.[71] obtained by numerical methods based upon SUSY. ( Note 
that the ISPP relations are preserved for arbitray value of £ /3' due to the 
'scaling' property given by eqn.(191) of the text ). 



]{AHO) E iDW 0) EiAHO) E (ow 0) 

ref.[71] ref.[71] 






1.95608 


2.38721 


1.93548 


1.93548 


1 


6.37732 


6.24897 


6.29849 


6.29849 


2 


11.7352 


11.3668 


11.6810 


11.6810 


3 


17.9931 


17.4785 


18.0426 


18.0426 


4 


25.0597 


24.4375 


25.2546 


25.2546 


5 


32.8581 


32.1484 


33.2261 


33.2261 


6 


41.3276 


40.5427 


41.8910 


41.8910 


7 


50.4197 


49.5679 


51.1979 


51.1979 


8 


60.0950 


59.1822 


61.1053 


61.1053 


9 


70.3204 


69.3513 


71.5790 


71.5790 


10 


81.0680 


80.0462 


82.5899 


82.5899 


11 


92.3136 


91.2421 


94.1129 


94.1129 


12 


104.036 


102.917 


106.126 


106.126 


13 


116.217 


115.053 


118.611 


118.611 


14 


128.839 


127.632 


131.549 


131.549 


15 


141.889 


140.640 


144.927 


144.927 


16 


155.351 


154.062 


158.728 


158.728 


17 


169.214 


167.887 


172.942 


172.942 


18 


183.467 


182.102 


187.557 


187.557 


19 


198.099 


196.698 


202.561 


202.561 
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7. Comparision of the Results of NGAS with 
the Exact Predictions of Super Symmetry for 
the Sextic Oscillator 



7.1 Super Symmetric Quantum Mechanics (SUSYQM) 
in Brief 

Supersymmetry (SUSY) is a symmetry beween Fermions and Bosons 6 
which was invoked by physicists to obtain a unified description of all basic 
interaction of nature although there has been no experimental evidence of 
SUSY being realised in nature. Nevertheless, in the last fifteen years, the 
ideas of supersymmetry have stimulated new approaches and applications 
to the different branches of physics including atomic, molecular, statistical 
and condensed matter physics as well as non-relativistic quantum mechanics 6 . 

In this section, we will be concerned with the application of SUSY in 
the area of non-relativistic quantum mechanics - a field now known as su- 
per symmetric quantum mechanics (SUSYQM). The investigations in this 
area have rapidly proliferated in the past two-decades which have estab- 
lished SUSYQM as an important and ineresting area of reaserch. The 
application of SUSY to quantum mehanics has yielded several new and sig- 
nificant results. In particular, the class of potentials yielding exact analytic 
solution of the non- relativistic Schrodinger's equation have been classified 
and understood readily on the basis of SUSY. Moeover, the results and con- 
cept of "partner-potentials" yielding iso-spctral hamiltonians are among the 
significant predictions of SUSYQM. For reviews of the subject, the ref.[6] 
may be consulted. 

Our immediate interest in the subject of SUSYQM arises in this thesis 
from the fact that the latter has exact-predictions of iso-spectrality of the 
sextic-AHO and the sextic-DWO which form partner hamiltonians for spe- 
cific choices of the parameters defining the potentials. To this topic, we turn 
to in the next section. 

7.2 Super Symmetry and Iso-spectrality of Partner 
Potentials 
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One of the simplest non-trivial applications of SUSYQM ' is made for 
the case of the sextic-oscillators (AHO and DWO). Consider the "super po- 
tential" ref.[6]: 

W(4>) = (3<t>\ (178) 

This is the simplest anharmonic, supersymmetric and "parity even" po- 
tential with no known analytic solutions. This is perhaps the simplest choice 
beyond the linear solvable case. Acording to SUSY quantum mechanics 6 the 
" partner-potentials" are given by: 

V { - ] = -{W 2 -W) (179) 

V { - +) = \{W 2 + W) (180) 

By using eqn.(178) in eqs.(179) and (180) the "partner-potentials" are 
generated as given below: 

v^ = l -{p<f-m 2 ) (isi) 

= h(3 2 ( f) (i + 3(3 ( j) 2 ) (182) 

Both the "partner-potentials" correspond to parity- even, confining po- 
tentials. One of these, corresponds to the double-well. The ground 
state energy of the "partner-potential" occurs at E = 0. In this case 
there is no trace of the double well structure as the zero-point energy is just 
equal to the well depth. 

The potentials for the sextic DWO and the AHO which are used in the 
Hamiltonian given by the eqn.(163) and eqn.(142) respectively can be written 
as: 

V 1>2 = A0 6 T ^0 2 (183) 

where A > and g > . On the other hand, the Hamiltonians 
corresponding to the "partner-potentials" are: 
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= \p 2 + 



(184) 



= -p 2 + (185) 
2 

In standard notation, the 'exact' results of SUSY for the above Hamilto- 
nians as given in eqs. (184-185) can be summarized as follows 6 : 

ii) Eit=Ei + \ (186) 

{ii) E^ = 0, (187) 

(in) = A e M- f W(y)dy), (188) 

where, n=0,l,2, ; the ground state wave function for if is denoted by 

ipQ \<p) and 'A' denotes its normalisation. The property given by eqn.(186) 
is referred to [6] as "Iso-spectrality" of Partner Potentials (ISPP). Eqn.(187) 
is a rigorous result of exact (unbroken) super symmetry, while eqn.(188) is 
the prediction for the ground state wave function of H^~\ 

Application of eqs. (181,182,184-188) to the case of the sextic AHO and 
DWO characterised by potentials in eqn.(183) becomes at once obvious, when 
the following specific values for 'A' and 'g' are chosen: 

A = g = 3(3; (3 > 0. (189) 

For the above choice of 'A' and l g\ eqn.(186) can then be rewritten as: 
E { n D + r 0) (X,g) = E^ \X,g) (190) 



7.3 Comparision of the Energy Levels obtained in 
NGAS (LO) with the Exact Results from SUSY 
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Surprisingly, it is found that the relation given in eqn.(186)(or eqn.(190)) 
is obeyed to a very good accuracy by the LO results from NGAS for all allowed 
values of A > 0! In Table-5, we demonstrate the (approximate) validity of 
the ISPP relation of SUSY in NGAS, by comparing the energy level of the 
DWO for the excitation label '(n+1)' with that of the AHO for the label 'n'. 
The agreement is seen to be impressive, particularly at large values of l n', 
considering that only the LO-results are used. It may be observed, in this 
context, that the formulae for energy levels of sextic oscillators in NGAS, 
given by eqs.(162) and (177), obey the following interesting "scaling" law: 



This scaling property guarantees the validity of the ISPP relation, eqn. (190) 
for arbitrary values of j3 , once the relation is established for any particular 
given value of the latter. 



7.4 Positivity Property of the Energy Eigenvalues of 
the Sextic-DWO predicted by SUSYQM and NGAS 



The other observation is regarding the "positivity" property of the energy- 
eigen values of the sextic-DWO predicted by SUSY through the eqs. (186- 
187), which is otherwise not obvious owing to the double- well structure of 
the potential (at least, this is not the case for the SSB-phase of the QDWO!). 
Interestingly, the positivity of the energy levels of the sextic-DWO, as pre- 
dicted by SUSY, is also dynamically realized in NGAS. This is because of 
the fact that the SSB-phase is ruled out on grounds of stability (see, remarks 
following eqn. (173)). 

As a final confirmation of consistency with the exact results from SUSY 
on energy levels of above systems, it is necessary to establish not only the 
ISPP -relation, but also the absolute magnitudes of the former. In Table- 
4(page-66), we compare the results of ref.[71], on the energy levels of sextic- 
oscillators obtained by sophisticated numerical methods, with those based 
upon the simple formulae, eqs. (162,177) in LO of NGAS for A = 0.5 and 
n < 20 . It can be seen from this comparison that there is good agreement 
for all values of 'n', except for the ground state-energy of the DWO. In the 
latter case there is some deviation from the exact result of SUSY, given by 
eqn. (146). It may be plausible that the discrepancy could be due to the 




(191) 
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departure from the predicted exact ground state wave function as given by 
eqn.(147) from the wave function of the DWO in LO of NGAS. This aspect 
is further investigated below: 

7.5 The DWO-Ground State Wave Function in NGAS 
and SUSYQM 

Having demonstrated the approximate validity in NGAS, of the ISPP re- 
lations and positivity property of energy levels predicted by SUSY, it remains 
to compare the respective ground state wave functions. The exact result from 
SUSY is given by eqn.(188). For the case of sextic DWO, this result can be 
made more specific as given below. Using eqn.(178) in eqn.(188) the ground 
state wave function is written as: 

4 _) = A expi-fi^/A) (192) 

The normalisation constant 'A' can be calculated easily and is given by 

A = (SflW/y/ni/A), (193) 

where T(z) is the Euler Gamma-function. Hence, the SUSY prediction for 
the ground state wave function is given by 

4 DWO \<I>^)\susy = m 1/8 exp(-^/4)/Vf(l/4) (194) 

Using the value of T(l/4) , eqn.(7.17) can be approximately represented 

as: 

4 DW °\<f>,f3)\susY = (O.681O8)/5 1 / 8 exp(-/30 4 /4). (195) 

On the other hand, the ground state wave function in LO-NGAS corresponds 
to that of an effective simple harmonic oscillator with variable frequency 
determined by the corresponding gap- equation. For the case of the sextic- 
DWO, the NGAS result is given by 

4 WO) (0^)k G A5 = (^) 1/4 e-^ 2 / 2 (196) 

7T 

where 'oV is the solution for the ground state of "gap-equation" given in 
eqn.(174) i.e., 

^ 4 + ^a-(y)A = (197) 
70 



The physical, acceptable solution of the above "gap-equation" is given by 

^[(vV + gOA)-^] 1 / 2 (198) 

On substitution of A = \0 2 \ g = 3(3 , the eqn. (198) leads the following 
result 

( ^0 )V 4 = (1^111)1/4^1/8 = . 828(3 l/8 (199) 
7T 71 

Using eqn. (7.22) in eqn. (7. 19) we have the following equation 

= 0.828f3 1/8 exp(-u a (j) 2 /2). (200) 

where 'u; a ' satisfies the gap-equation, eqn. (197); is defined by eqn. (189) 
and £ = 1/2 , corresponding to the ground state of the DWO. Comparision 
of the coefficients of the exponential terms in eqs. (194) and (200) interest- 
ingly reveals that not only the "/3 1 / 8 " factor is common but also that the 
coefficients: 0.68108 and 0.828 are comparable. We compare the two results 
in Figure 1 for (3 = 100. The quality of the approximation can be judged 
from this figure. It is plausible that the inclusion of higher order corrections 
to the ground state wave function in IPT of NGAS (Chapter-10) may fur- 
ther improve the agreement. 

To summarize the results of this Chapter, it is shown that NGAS re- 
spects and preserves the exact results of SUSYQM with good accuracy. 
It would be interesting to extend the comparison 120 to the system of self- 
interactin oscillators described in SUSYQM by the super potentials: 
W± = [3(fi 3 ± 70 , which generate a family of sextic-, quartic- and quadratic- 
AHO/DWO for different values of the parameters '/?' and '7'. This is, how- 
ever, beyond the scope of the present thesis. 
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i HtW- 1 1 r 

- 100 -0-60 -0,20 030 0.60 1.00 

Figure 1: 

Comparision of the ground state wave function of the sextic-DWO predicted 
by SUSY ( curve with sharper peak ) with that obtained in LO of NGAS for 
P = 100, see eqs.(195) and (200) of text. 
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8. The NGAS for the Octic-Anharmonic 
Oscillator ( Leading Order(LO) Results ) 



8.1 The Octic-Anharmonic Oscillator: 
Introduction 

The octic-anharmonic oscillator, like its quartic- and sextic- counter- 
parts, finds applications in modelling molecular physics, lattice-vibrations 
in solids and in quantum chemistry. Because of the higher anharmonicity, 
the system also provides the theoretical laboratory for more stringent-tests 
for non-perturbative approximation schemes in quantum theory since the di- 
vergence of the naive ( Rayleigh-Schrodinger) perturbation theory becomes 
still more severe 12 ' 13 ' 14 in this case. 

The system is considered here to test the generality of the application of 
NGAS to the case of still higher anharmonicity and hence to test the relia- 
bility of the said scheme. The detailed application is described below. 

To demonstrate further the generality and uniformity of the approxi- 
mation (NGAS), we apply the method to the case of the next higher an- 
harmonicity, i.e., the octic-anharmonic oscillator, described by the following 
Hamiltonian: 

8.2 Application of NGAS to the Octic-AHO 

8.2.1 The derivation of the Gap-equation for the Octic-AHO 

As before, we start from the Hamiltonian describing the system: 

H =\v 2 + \gtf + A0 8 ; g, A > 0. (201) 

The free field Hamiltonian corresponds to H s = |p 2 + \g<fi 2 and the 
interaction term is given by AHj((f)) = A0 8 . To develop the NGAS for the 
octic AHO we follow the identical ansatz for V(<f>) as given in eqn.(62): 

V(<P) = A<f - B<f) + C (202) 
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The parameters A,B,C can be determined self-consistently as before. For 
this purpose, consider the "Effective Hamiltonian (EH)" in this case given 
by, ' 

#oOM = \p 2 + \g<P 2 + w{<p) (203) 

As in the case of quartic- and sextic- anharmonicity, eqn.(203) is trans- 
formed into the following diagonal structure by using the eqn. (202) and given 
by: 

#o(<M = l -p 2 + \u 2 {<P-a) 2 + h Q (204) 

where again, 

u 2 = g + 2\A, (205) 



a = \B/u 2 , (206) 

ho = XC - \uj 2 a 2 . (207) 

Eqn. (204) corresponds to "shifted" effective Harmonic oscillator. The 
field '0' is shifted by 'a' and the energy is shifted by l h \ As in other cases 
here also the physical requirement is that uj > ; a = real. 



Using the creation- and annihilation operators defined by the eqs.(67) and 
(68) along with the equal-time commutation relation defined by the eqn. (69) 
the EH, H given in eqn. (204) is expressed into the desired diagonal form 

Ho = uj (N b + 1/2) + h . (208) 

To obtain eqn. (208), we have introduced the number operator N b = tfb 
and its eigen-states by, N b \n > = n\n >, < m\n > = 5 mn as usual. 
This leads to the energy spectrum as in the earlier cases: 

= u^ + ho, (209) 

where f = (n + 1/2) ; n = 0,1,2, 
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It is next required to determine the frequency l oo 1 and '/iq' defined by 
the eqn.(207). For this purpose, we note that the quantum average of the 
eqn.(201) is given by: 

<H> = <H > = \<P 2 > +\g < 2 > +A < 8 > (210) 

where eqn.(55) has been used and < >, <0 2 >, < p 2 > 
have been calculated using the standard properties of creation- / annihilation 
operators which are given in eqs.(80) and (81). The QA of the last term in 
eqn.(210) can be obtained by standard procedure and is given by: 

< 8 >= a 8 + 28a 6 (i) + ^ 4 (4£ 2 + 1) + f (%4£ 2 + 5) + 

^(-)(16£ 4 + 56£ 2 + 9) (211) 

Substituting the QA- values in eqn.(210) we obtain the equation 

<H > ~y + + A[a 8 + 28a 6 (i) + ^<r 4 (4£ 2 + 1) + 

2 2 UJ 4u; z 

f (^)(4^ 2 + 5) + |^)(16£ 4 + 56C 2 + 9)] (212) 

Applying the variational minimisation condition d < H > /duj — 
leads to the following equation 

uj 5 ~uj 3 (g + 5QXa 6 )-105uj 2 (Xa 4 /^(4^ + l)-105ujXa 2 (4^ + 5)-35Xh^) = 0, 

(213) 

where h(£) = £ 3 + (7f/2) + (9/160- Eqn.(213) is referred as the "gap- 
equation" for the case of octic AHO. 

Again carrying out minimisation with respect to a i.e., d < H > /da = 
we have the following equation 



a [g + X(8a 6 + lQ8a A ^/uj) + 105a 2 (A^ 2 + l)/uj 2 + 35^ 2 + 5)/uj 3 ) ] = 0. (214) 

Eqn.(214) corresponds to equation to ground state (EGS) for this case. 
The solution of eqn.(214) for a leads the two following equations: 

(i) a = 0. (215) 
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(ii) [^ + A(8 ( t 6 + 168ct 4 (^) + 105 ( t 2 (4^ 2 + 1)/cj 2 + 35^(4^ 2 + 5)/cj 3 ) ] = 0. 

(216) 

Obviously, eqn. (216) has no physically acceptable solution for A, g, uj > 0. 
Hence the 'physical' solution of the EGS, eqn. (214), is at a = 0. Substitu- 
tion of this value in eqn. (213) leads to the simplified GE, given by: 

uj 5 - gto 3 - 35A/i(0 = (217) 

It is then straight forward to evaluate A, B, C appearing in the approx- 
imating potential, eqn. (202) in terms of 'u;' and l a\ These are given by the 
following expressions: 

fi 105a 4 (l + 4£ 2 ) ,105a 2 ^ ,_ 2N 35/i(£) , rt , n , 

A = 28ff + — kr^ +( ^ )(5+4a+ ^' (218) 



B = a[(l+9) i v+8wV 6 + 168<rV+105 t T 2 (l+4f 2 ) + (^)(5+4^)], (219) 



C = <<p 8 >-A<(j) 2 >+B<(j)>. (220) 



8.2.2 Solution of the Gap Equation and Determination of the 
Energy Spectrum 

The solution of eqn. (217) determinines the frequency 'cj' of the "shifted" 
harmonic oscillator. To obtain the energy levels one substitutes a = and 
'uj^ as the solution of eqn. (217). This leads, after some simplification, to the 
following simple formula: 

E^\ ocUc -aho = d)(5cu + ^), g > 0. (221) 

o UJ 

where 'a;' is obtained by solving eqn. (217) numerically. 

In Table-5 ,we compare this LO- result in NGAS with earlier computations 133 
over a wide range of values of 'A' and 'n'. It can be seen from this comparision 
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that the results obtained in the LO of NGAS are already quite accurate over 
the full range of the parameters, which demonstrates the generality of the 
method and uniformity of the approximation with increasing anharmonicity. 

We next turn our attention to the physics of the effective vacuum state: 
\vac > , obtained as an approximation to the true vacuum of the theory. 
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Table- 5: Sample results for the octic- AHO in the LO NGAS compared 
with results of earlier calculations from ref.[133] (shown in parentheses) , over 
a wide range of '\ ' and 'n'. 



n 


A = 0.1 


1.0 


5.0 


50.0 


200.0 





1.3005 


1.7794 


2.3290 


3.5565 


4.6425 




(1.2410) 


(1.6413) 


(2.1145) 


(3.1886) 


(4.1461) 


1 


4.4717 


6.3946 


8.5167 


13.172 


17.259 




(4.2754) 


(5.9996) 


(7.9296) 


(12.1950) 


(15.9519) 


2 


8.6264 


12.717 


17.126 


26.698 


35.062 




(8.4530) 


(12.421) 


(16.711) 


(26.033) 


(34.183) 


4 


19.763 


30.026 


40.863 


64.165 


84.444 




(19.9930) 


(30.4605) 


(41.495) 


(65.202) 


(85.8251) 


6 


34.217 


52.669 


72.044 


113.48 


149.47 




(35.0560) 


(54.1403) 


(74.0830) 


(116.7629) 


(153.83) 


8 


51.570 


80.013 


109.65 


172.99 


227.97 




(53.146) 


(82.6496) 


(113.3486) 


(178.9215) 


(235.82) 


9 


61.239 


95.255 


130.64 


206.23 


271.81 




(63.225) 


(98.5529) 


(135.26) 


(213.6157) 


(281.5864) 


10 


71.532 


111.49 


153.01 


242.64 


318.52 




(73.954) 


(115.49) 


(158.5991) 


(250.5751) 


(330.3433) 


11 


82.424 


128.68 


176.69 


279.14 


368.06 




(85.308) 


(133.42) 


(183.3103) 


(289.71) 


(381.97) 


12 


93.893 


146.79 


201.65 


318.67 


420.14 




(97.2636) 


(152.31) 


(209.3443) 


(330.9440) 


(436.3695) 


13 


105.92 


165.79 


227.84 


360.14 


474.85 




(109.7967) 


(172.11) 


(236.6436) 


(374.1834) 


(493.4143) 


14 


118.49 


185.65 


78 

255.21 


403.50 


532.06 




(122.89) 


(192.81) 


(265.1732) 


(419.3737) 


(553.0335) 



9. The Property, Structure, Stability and 
the Significance of the " Effective" Vacuum in 
NGAS 

9.1 The Bogoliubov Transformation relating the Free- 
Field Vacuum to the Effective Vacuum 

The study of the properties and the structure of the vacuum of interacting 
quantum systems are of considerable importance 146 . In the present scheme, 
the vacuum state, \vac > of the effective Hamiltonian H , approximates 
the vacuum of the true interacting theory in the leading order. To study its 
properties and structure in comparision to the "free"- field vacuum, |0 > , 
it is useful to start with eqs. (71-73, 88-90). In view of eqn.(90), the creation- 
and annihilation operators of the 'free-theory' and the approximated the- 
ory with self-interaction, are related by quantum-canonical transformation 
("Bogoliubov- Transformation") 52 , given by: 

b = a cosh(a) — a) sinh(o;) (222) 

¥ = a) cosh(a) — a sinh(o;) (223) 
The two vacua are then related by the following equations: 

\v ac>= exp[(l/2) tanh(a) (aV - aa)} |0 > = U(a; a, a f )|0 > (224) 

The derivation of eqn.(224) follows from eqn.(222) by using the defining 
property of the vacuum, b\vac >= and the representation of the annihila- 
tion operator given by, a = 5 /5a) . The parameter 'a' occuring in the above 
equations, can be simply related to 'id' by using the eqs. (71,88,222) as given 
below: 



From eqs. (71) and (88) one can obtain, 

b + &t = J^- (a + a f ) (225) 
V 

Again from eqs. (222) and (223) we have the relation 
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b + &t = ( a + a t) e -a (226) 

By solving eqs.(225) and (226) we get the relation between 'a' and 'u;' 
and is given by 




(227) 



Then, 

a = (1/2) ln(u /u), uj = y/g (228) 

It is useful to have the transformation inverse to eqs. (222-223). This is 
given by 

a = b cosh(ct) + b ] sinh(a) (229) 

a) = &t C osh(a) + 6 sinh(a) (230) 

The following significant physical results follow from the above equations, 
eqs. (222-225). 

9.2 Structure of the effective vacuum 

(i) A non-trivial structure ("dressing") of the "effective" vacuum(EV) of 
the theory emerges from the equations. The situation could be analogous to 
the case of the ground state of the super-fluid 147 and the hard sphere Bose 
-gas 49 . The structure is characterized by the non-vanshing number density 
of the free particles in the EV and its critical dependence on the strength of 
the interaction, which is shown below: 

The number density of the free particles in the effective vacuum (EV) is 
defined by the relation 

n = < vac\a) a\vac > (231) 

From eqs. (229) and (230) and using the defining property of the vacuum, 
b\vac > = = < vac\$ one obtains the relation 

< vac\a) a\vac > = sinh 2 (a) (232) 
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Hence by using eqn.(227) we have the following equation 



n = < vac\a)a\v ac > = sinh 2 (a:) = -( 1 - — 2) (233) 

4 UJq uj 

Considering the case of the quartic- AHO, it can then be shown using the 
'gap- equation' that no ~ A 1 / 3 for A >> 1 . In the limit of vanishing 
interaction, one recovers the expected behaviour, n — > for A — > . 

(ii) Secondly, eqs. (222-224,228) imply an entirely new physical interpreta- 
tion of the parameter, ' oo ' which determines (through eqs. (228) and (235)) 
the "vacuum structure function" 'a' in the sense that q ^ (i.e., u ^ ) 
signifies the non-trivial structure of the EV in presence of interaction. 

9.3 Stability of the Effective Vacuum 

In the remaining part of this Chapter, we investigate the stability prop- 
erties of the EV. 



Instability of the Perturbative( "free-field") Vacuum 

It is shown below that the perturbative vacuum, |0 > becomes unstable 
compared to the effective vacuum \vac > for all values of the coupling 
strenth 'A'. For this demonstration we consider, for reasons of simplicity, the 
case of the QAHO. The standard method for studying the stability properties 
is to consider the "effective potential (EP)" . The EP, for any given choice of 
a vacuum state, is defined 148 to be the expectation value of the Hamiltonian 
in the chosen vacuum-state and expressed as a function of the VEV of the 
"field" '(f)'. For the case of the QAHO, this is obtained in LO of NGAS, 
from eqn.(86) by choosing, g — l,n — and 'a;' constrained to satisfy 
eqn.(92). The resulting expression is as follows: 



where, 



ngas / \ u , (1 + 12Aa 2 ) 3A 

eff ((T) = I + 4^ + 4^ + K ' (234) 



V c = l -a 2 + \a\ (235) 
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is the 11 classical potential' and 'u;' satisfies eqn.(92). 

An analogous expression for the corresponding EP based upon the perturba- 
tive (free- field) vacuum is obtained by the substitution, uj — > 1 in eqn. (234) 
above ( this follows by comparing, eqs.(71) and (88)) and is given by 

V e p f f(a) = \ + 3A(a 2 + \)+V e (236) 

The ground state energy is defined to be the global minimum of the effective 
potential and corresponds to a = in either case. We thus obtain 
the respective ground state energies from eqn. (234) given by the following 
equations : 

E » = i + i + i7 = (1/8)(3 " +1/w) (237) 

The above relation is obtained by using the GE for the ground state: 
uj 3 — uj — 6A = . However, E Pert is obtained from eqn. (236), by using 
o" = and given by 

E Pert = 1/2 + 3A/4. (238) 

(Note that eqn. (237) is also contained in eqn. (109) for the special case 
considered here (i.e. g — 1, n = )). Recalling that the GE for the ground 
state is given by: uj 3 — uj — 6A = , it is straight forward to establish that: 

E - E Pert < 0, for all values of A. (239) 
which is shown below. 

By using eqs.(237) and (238) we have the relation 

1 1 1 

AE(X) =E - E^ = (a; + -) + -(-- 1) - 1/2. (240) 

4 UJ 4 UJ 

Here it is to be noted that as A — > ; uj — > 1 then AE(\) — > as 
it should! Now let us study how uj(\) and AE(X) behave as A ~ + . 
Recalling the gap equation for the ground state we have the equation : 

uj 3 - uj - 6A = (241) 
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By substituting to = 1 + £(A) ; |£(A)| << 1 in eqn.(241) then 
£ ~ 3A + 0(A 2 ) as A -> . Hence in the lim A ^ , u(X) ~ 1 + 3A. 
So it is seen that 'a;' is monotonically increasing with 'A' as 'A' increases! 
Substituting this value of 'a;' in the eqn.(240) and after some simplification 
we get the relation 

q\2 

lim AE(\) = E - = -(— ) < (242) 

From this relation it is evident that AE(X) is monotonically decreasing 
with increasing 'A' (for A > ) ! This proves the instability of the pertur- 
bative (free-field) vacuum of the QAHO. It may be noted that, although we 
have established the above result for the QAHO, the same can be rigorously 
demonstrated in all other cases of anharmonicity considered here. 

9.4 Significance of the Effective Vacuum state 

In the context of quantum-mechanics, the results presented in the previ- 
ous sections: 9.1 - 9.3 are significant in many respects which are listed below: 

(a) In the first-place, the results contained in eqs. (231-233) demonstrate that 
the effective vacuum state is endowed with non-trivial particle-content and 
structure, which is directly attributable to interaction since the particle- 
content vanishes as A — > . This result is analogous to the "dressing" 
of the "physical" -vacuum and is made more transparent when one considers 
the A0 4 -field theory discussed subsequently in Chapter-11. 

(b) The other significant result is the equivalence of the NGAS- vacuum with 
that obtained through quantum canonical/ Bogoliubov-transformation as de- 
picted in eqs. (222-230). It also demonstrates that the effective- vacuum 
(EV) state is not obtainable in perturbation-theory due to non- analytic 
dependence of the "vacuum-structure function (VSF)" 'a' on the coupling- 
strength £ A'( see,eqn.(228) and the gap-equation, e.g., eqn.(241)). 

(c) The third significant observation is contained in the conclusive demon- 
stration of the instability of the "perturbative" / free-field vacuum in respect 
of the EV, as contained in eqn.( 242). In this context it may be further sig- 
nificant that the convergence of the perturbation theory about the free-field 
vacuum discussed in the next Chapter may have a direct bearing with the 
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stability properties of the theory. This is discussed in greater-detail in the 
next Chapter. 



10. Improved Perturbation Theory(IPT) in 
NGAS 

10.1 Considerations of the Convergence of the 'naive' 
Perturbation Theory and other Variants 

One of the main motivations for proposing the NGAS as described earlier, 
is the possibility of construction of an improved perturbation theory (IPT) 
which could be convergent for all allowed values of of the coupling strenth, 
'g' and 'A'. This expectation is based upon the result, eqn.(59), which is 
reproduced below: 

< n\XH'\n >= 0. (243) 

Since H = H + XH' , eqn.(243) naturally suggests that the IPT 
be constructed by choosing Hq as the unperturbed Hamiltonian and XH' 
as the perturbation. The convergence of the resulting IPT is intuitively 
suggested since the basic condition of convergence 149 is satisfied by ensuring 
that the magnitude of the perturbation remains always sub-dominant to the 
unperturbed contribution: 

| < n\XH'\n > | = << | < n\H \n > | (244) 

The important point to note is that eqn. ( 244 ) holds for arbitrary values 
of 'q\ 'A' and V. In this context, it may be noted that the analogous re- 
quirement, which is the necessary condition for convergence of perturbation 
expansion, does not hold good in the case of naive perturbation theory 
(NPT), where the entire self- interaction, XHj((p) is chosen as the perturba- 
tion to the '/ree' Hamiltonian, H s (p,<f)) (see eqn. (52)). Consequently, the 
divergence 12,13,87 of the NPT is anticipated as it could be traced 14 to the 
eventual dominance of the perturbation-contribution over the unperturbed 
one for any value of A > 0, no matter however small. This was explicitly 
demonstrated in ref. [12,13,87]. 
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In the next section, we formulate the improvement of the LO results 
through the "( improved ) perturbation theory" (IPT) developed using XH' 
as the perturbation while treating H as the unperturbed Hamiltonian. 

10.2 Improvement of the LO Results through the IPT 
in NGAS 

The unique feature of NGAS summarized in eqn.(244) leads to the sys- 
tematic further (order-by-order) improvement of the LO results (which are 
already accurate to within a few percent of the exact result). We demon- 
strate below, the improvement in accuracy, by inclusion of the higher-order 
contribution in IPT for the case of the QAHO and the QDWO. 

In the Rayleigh-Schrodinger (RS) development of the perturbation series, 
the perturbative correction to the energy levels is given by the standard 
expansion: 

E n = £f + A£« + + + (245) 

where, the LO-contribution, E^> has already been defined, ( see, eqs.(54- 
55)). It is important to note that the first order contribution AE^ vanishes 
due to eqn.(243): 

AE^ =< n\XH'\n >= 0, (246) 



(In the above sense, the IPT can be regarded as optimal and this result, 
eqn.(246), distinguishes the IPT from many other variants of perturbation 
theory 150 used earlier, for the problem). Using eqn.(246), the next higher 
order (HO) contributions are given by the following expressions: 

A£f = £ l(A#;jl7A„ m , (247) 

AE^ = £ (^»m) i^H'mk) i^kn) ^48) 

Similar expressions for still higher-order corrections can be obtained by standard 149 
methods. In the above equations, we have used the following notations: 
(XH') mn =< m\XH'\n > and A nm = (E$ - E^). The above formulae can 



85 



be applied to the case of the QAHO as discussed below. 

10.3 IPT- Applied to Quartic-Oscillator (AHO/DWO) 

The Hamiltonian for the systems given in eqn.(61) is 

H = \p 2 + ^0 2 + A0 4 , (249) 
The free field Hamiltonian is given by 

HM) = \p 2 + I 9 4? (250) 
and the self-interaction term is given by 

A# 7 (0) = A0 4 . (251) 
With the new ansatz V(4>) the effective Hamiltonian is given by 

#oOM = \v 2 + \g<P 2 + W{<P) (252) 
where V{4>) is given by the equation 

V ((f)) = A <f) 2 - B <f) + C. (253) 

The parameters A,B,C are determined and given by the ens. (95), (97), 
(99) with the choice a = (for the case of the QAHO). The modified 
interaction XH' = H — H = (A0 4 — V(<f))) have vanishing quantum 
average for arbitrary 'A ', '<?' and 'n' as has been noted. Hence, the first order 
contribution vanishes as given in eqn.(246). To calculate the higher order 
perturbative contribution, the matrix elements for the case of the QAHO, 
are given by 

(XH') mn =< m\X(t) A \n > -(3A/w)/(f) < m|0 2 |n >, m ^ n (254) 
where, 

< m\4> 2 \n >= (l/2u) (5 m , n +2 ^(n + 1 ) (n + 2) + 

$m,n-2\fn(n - 1)). (255) 
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and 



< m|0 4 |n >= {l/Auo 2 ) {5 rn . n+ ^ {n + 1) (n + 2) (n + 3) (n + 4) + 

S m , n -4\Jn(n -l){n- 2)(n - 3) + 
(2n + 3)^/(n + l)(n + 2)) + 
2(5 m ,„_ 2 (2n - l)^n(n-l))) (256) 

(The above matrix elements are easily calculated by using eqn.(71) and 
introducing the basis states \vac > , and \n > ; n = 1,2,3.. with the 
defining property of the effective vacuum b\vac > = = < vac\b\ and 
further relations: \n > = \vac > ; N = tfb ; N\n > = n\n > ; 

tf\n> = \/n + l\n > and b\n > = ^pa\n — 1 > ). 

We present in Table-l,(see page -46) the results for the energy levels of 
the QAHO, with the inclusion of the second-order perturbation correction. 
In the same Table we also compare our results with available 'exact' numer- 
ical results and results of calculation in second order perturbation theory of 
ref.[44], which is based upon the operator methods. It may be seen from 
this Table that the accuracy is considerably improved by inclusion of the 
perturbation correction and further that the convergence of the IPT is found 
superior (order- by-order) to that in ref.[44]. Similar results for the QDWO 
after inclusion of the second order perturbative correction in IPT, is pre- 
sented in Table-2 (at page-54). In this Table, we also compare our results 
with those obtained by inclusion of twenty orders in the 'modified' pertur- 
bation theory of ref. [36]. Again uniform improvement in accuracy is seen. 
It is also seen from the Table that the results of the present analysis which 
includes only the first non-trivial perturbative correction, compares well with 
the results of ref. [36] obtained by sophisticated numerical methods. In the 
context of the above results, the following observations may be relevant: 

10.4 Remarks on Convergence and other aspects of 
IPT in NGAS 

(a) Corrections up to the fourth-order in IPT have been computed for 
the QAHO and the QDWO although we have reported only the second order 
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correction in the Tables-1,2. It is seen that these higher order corrections 
remain uniformly small compared to the LO results over the full range of '\ ' 
and e n' and decrease fast with the order of correction, which is consistent with 
the expectations from a rapidly converging sequence. 

(b) As has been demonstrated in the previous section, the 'perturbative' 
ground state (i.e., corresponding to the free-field Hamiltonian H s (<f),p) , see 
eqn.(52)) becomes unstable compared to the ground state of the 'effective' 
Hamiltonian, H . Thus the stability of the theory, and the convergence 
of the IPT- both appear to critically depend on the choice of H . It may 
perhaps be plausible, therefore, to conjecture that the convergence of the 
perturbation theory may be intimately connected with the choice of a stable 
vacuum resulting from a proper effective Hamiltonian chosen as the unper- 
turbed part. 

(c) When compared with the results of some other variants of perturba- 
tion theories 150 applied to the above systems of anharmonic and double well 
oscillators, the IPT appears to provide better convergence, when compared 
at each order. 

(d) For the case of the QDWO, the results of IPT, as well as those from 
the other variants of perturbation theories, show poor convergence near the 
transition point: A ~ A c (£) as intuitively expected 151 . However, since A c is 
small (A c (£) < 0.0362886) , this limitation does not affect most applictions 
of practical interest. In particular, the strong coupling regime, A >> 1 is 
excellently described by the IPT. 

(e) Although we have provided only plausible argument in support of the con- 
vergence of the IPT, a formal proof can be attempted following the methods 
available in the literature 152 . 

In the following Chapter, we extend the method to A0 4 - quantum field 
theory in ( 3 + 1 ) dimensions. 

11. Application of NGAS to A</> 4 - Field 
Theory in ( 3 + 1 ) dimensions 
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11.1 Formulation of the NGAS for A0 4 quantum field 
theory 

We consider, in this chapter, the application of the NGAS to A0 4 the- 
ory in (3+1) dimensions. In view of the successful application of the NGAS 
to the AHO's and DWO's as described in earlier chapters, it is but natural 
to extend the formalism to A0 4 quantum field-theory and test the conse- 
quences of the scheme. 

It may be worth while pointing out here that the A0 4 - field theory in 
physical dimensions ( and in lower dimensions ) is an important physical sys- 
tem which finds crucial applications in diverse areas of physics, e.g., the stan- 
dard model of particle-physics 84 , cosmology 153 , condensed matter physics 121 , 
phase-transitions and critical phenomena 75 etc. Besides, this theory pro- 
vides the simplest theoretical laboratory for testing the various approxima- 
tion schemes in quantum field-theory(QFT). Hence it becomes imperative to 
test the current approximation scheme NGAS by applying the same to A0 4 - 
QFT. 

In this thesis, we consider the theory in the massive, symmetric-phase 
described by the Lagrangian: 

C = \{d,<t>){^4>) - \™ 2 <t> 2 - A0 4 , (257) 

where m 2 > . The Hamiltonian density derived from the above 
Lagrangian is given by 

H = i(m 2 2 + 2 + 2 ) + A0 4 (258) 
where we have defined: 4> t = d(f)(x,t)/dt and 4> a = d(j)(x ) t)/dx a . 

To formulate the NGAS for the above theory we follow analogous steps as 
in the cases of AHO/DWO considered earlier and choose an approximating 
potential (AP) denoted by V(<f)) such that ideally the defining constraints 
as given by eqs. (52-57) are satisfied. However, unlike the case in quantum 



89 



mechanics, the PEQA involving multiparticle-states, is hard to implement in 
QFT. Hence, on grounds of simplicity, we relax the condition, eqn.(56) by 
restricting the QA to be evaluated in the " few-particle" states only. To be 
more specific, we proceed as follows : 

11.1.1 Choice of V(</>) 

Using the most general ansatz quadratic in the fields 0, we parametrise 
V((f>) as given below : 

V(<j>) = A(f) 2 - Bcf) + C (259) 
As before, the " effective Hamiltonian" (EP) H is defined as : 

Ho = i(m 2 2 + 0? + 0D + A\/(0) (260) 

Substitution of eqn.(259) in eqn.(260) leads, after some simplification, 
to the following expression : 

m = \M 2 e + \e a + (26i) 

where, 

Z(x,t) = <j>{$,t)-a , (262) 




(263) 



M 2 = m 2 + 2\A , (264) 

and 

h = \C-\m 2 cj 2 . (265) 

In addition, £ t = d£/dt ; £ a = d^/dx a etc. Eqn.(261) is atonce identified to 
be the Hamiltonian density of the hermitian scalar-field £(x, t) . This is not 
surprising since the AP was chosen accordingly. However, the important 
point to emphasize is that the AP is proposed to incorporate the effects of 
self-interaction even though the EH corresponds to that of an exactly solv- 
able system. To demonstrate this aspect it is first necessary to obtain the 
spectrum of H which is done as follows. 
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The diagonalisation of the EH given by eqn.(261) is straightforward by 
using the Fourier expansion in terms of creation- and annihilation operators 

£(x,t) = <f>(x,t)-a = /^[^V- + ^)e-], (266) 
where 

Vl k (M) = 2(27r) 3 v /|i 2 | +M 2 = 2 (2tt) 3 u k (M) , (267) 

and kx = k°t — k.x, as usual. The operators b(k), tf(h) satisfy the standard 
(equal-time) commutation relations (ETCR) : 

[b(k),b\$] = n k (M)5 3 (k-q), (268) 

which is a consequence of the ETCR between the 'field' <p(x, t) and its 
canonical conjugate momentum : ir(x,t) = dC/d(f> , given by : 

[<f>(x,t), n(y,t)\ = t5 3 (x-y). (269) 

11.1.2 Diagonalisation of the Effective Hamiltonian 

The energy of the system described by H is obtained by standard 
methods and given by : 

H = J d 3 xH (x,t) = l -J b\k)+b\k) b(k)}+ J d 3 xh 

(270) 

The spectrum of the states are analogously obtained and denoted by : | vac > 

I P >> I V\i P~2 >, etc where the effective vacuum state | vac > is 

defined by 

b(k) \vac > = , (271) 

and the multi particle-states are generated by multiple application of the 
creation-operator b\p) on | vac > : 

b\p) \vac > = | p > (272) 
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- = I Pi, P2 >, etc (273) 



V2! 

Let us note that, H Q \ vac > — E | vac > ; H \ p > = Ei (p) \ p > 
.....etc. where E , E\ etc correspond to the energy of the corresponding 
states. 

The next step is the implementation of the PEQA. 

11.1.3 PEQA for A0 4 QFT in NGAS 

This requirement translates to the following constraints : 

< vac\ 4 \vac > = < vac\ V(4>) \vac > , (274) 



<p\<j) 4 \p> = <p\V((f))\p>, (275) 

and similarly for multi-particle states. The implementation of eqs.(274, 275) 
require the evaluation of the QA of monomials of the field <p(x, t) such as : 
< vac\ (f) n (x,t) \vac > = < (f) n (x,t) >, < p\ 4> n (x,t) \p > etc. 
We first turn to evaluation of < (f) n (x,t) > . This is readly done using 
translational — invariance of the vacuum-state : \vac > and the ETCR, as 
given by eqn.(268). Some useful results thus obtained are given below : 

< <j>(x,t) > = a , (276) 



=° 2 + Jim ^ 2 + '°- < 277 > 

< 4 (f , t)> = a 4 + 6a 2 1 + 3/ 2 , etc. (278) 

Similarly, 

= ifSf) m (279) 

< ^ t)> = ItSh)^ m) (280) 



92 



11.1.4 Determination of the AP 

Using the above results in eqn.(274) one determines the AP uniquely, 
which is given by 

V{<t>) = 6(a 2 + J o )0 2 - 2a 3 (f> - 3(a 2 + I ) 2 . (281) 

In other words, the coefficients A, B, C defining the AP (see, eqn.(259)), are 
uniquely determined : 

A = 6(a 2 + J ) (282) 
B = 2a 3 (283) 

C = - 3(a 2 + J ) 2 . (284) 

It is further verified by explicit calculation that the above choice of the 
AP guarantees not only the equality of the QA over the vacuum state; i.e. 
< V{4>) >=< 4 > ( eqn. (274)), but also that for the one-particle states as 
well, i.e. < p \ V(4>) \ p > = < p \4> 4 \p >, as required by eqn. (275). 

It is important to emphasize at this point that eqns. (282-284) when con- 
sidered together with eqns. (263-265) form a complete set of self — consistency 
conditions which uniquely specify the physical consequences of the theory 
inthe leading — order (LO), such as the spectrum, renormalisation, stability 
properties and the structure of the effective vacuum. These physical conse- 
quences of the theory are discussed in the following sections. 

11.2 The "Effective Potential (EP)" and Renormalisa- 
tion in LO 

Before we discuss the ( non-perturbative) renormalisation programme in 
LO, it is useful to first investigate the consequences of eqs. (262-265) con- 
sidered together with eqs. (282-284). Substitution of eqn. (282) in eqn. (264) 
leads to the following equation : 

M 2 (A, a) = m 2 + 12Xa 2 + 12A/ (M 2 ) , (285) 
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i.e., 

M 2 (A, a) = m 2 + m^x 2 + 6A / f (fc) , (286) 

7 (27r) 3 v /|fc| 2 + M 2 

Eqn.(286) can be interpreted as the generation of the 'mass-gap' (i.e. shift 
in the bare-mass) due to interaction. In analogy with the terminology used 
earlier, we refer eqs. (285-286) as the "gap-equation(GE)" of the theory. This 
equation plays crucial-role in the subsequent discussions. 



Similarly, consideration of eqn.(283) together with eqn.(263) leads to the 
" equation for the ground-state (EGS)" 

°- [ <? " ^ ] = (287) 

As in case of the QAHO and QDWO's considered in earlier chapters, it 
is convenient to first obtain the solution of the EGS, eqn.(287). The two 
solutions of eqn.(287) are given by : 

(j = 0, (288) 

and 

^ ~ £ <-) 

We show below that eqn.(289) is not an acceptable solution on physical 
grounds. To establish this result substitute eqn.(289) in eqn. (285) which 
leads to : 

5M 2 = - (m 2 + 12A/ (M 2 )) (290) 

Since M 2 > is the physical requirement for definition of the theory 
( see, e.g. eqn.( 267 ) and eqn.( 270 )), equation (290) can not lead to 
acceptable solution unless 

A < i.e. A ee -g , g > ; (291) 

and 

12 gI (M 2 ) - m 2 > (292) 
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( Note that A < can not apriori be ruled out since bare-parameters 
in the Lagrangian are unobservable, see later). However, now the defining 
equation, eqn.(289) can then be rewritten as : 

a 2 = - — < , (293) 

which is not acceptable on grounds of the hermiticity i.e. 0t = 0! Thus 
eqn. (289) is ruled out as a solution and the unique physical solution of the 
ground state corresponds to eqn. (288). Having thus fixed the ground-state 
configuration, the implementation of the renormalisation programme can be 
done by defining the effective-potential (EP). The latter is defined as : 

U (a) = < vac\ H \vac > , (294) 

such that 

a = < vac\ (p(x,t) \vac > . (295) 

It may be noted that the l.h.s. of eqn. (294) is defined to be a function 
of a alone. This means that any other parameter occuring in 7io is to 
be variationally fixed by minimisation of < Ho >• The procedure is made 
explicit below by working out the current example, with Ho defined in 
eqn. (260). To this end, we first calculate < H > by using the eqn. (274) 
which guarantees the following equation: 

< Ho > = < H > = ]-m 2 < <f) 2 > + ]- < 4> 2 t > 

+ \<€> + A<0 4 >. (296) 

This works out to be : 



<n > = i raV + /o) + i/ i g y (^) + |^) 



m 

+ A ( <t 4 + 6a 2 I + 3I 2 ) . (297) 



This can be rewritten as, ( using ul(M 2 ) = \k\ 2 + M 2 ) 
< H > = h - \l ( m 2 + 12Acr 2 + 12A/ ) 

+ \m 2 ( a 2 + J ) + A ( a 4 + 6a 2 I + 3/ 2 ) , (298) 
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where, we have defined : 




= 0, ± 1, ± 2,.... (299) 



These integrals were first introduced by Stevenson 154 . Eqn.(298) can be 
simplified further by using the "gap-equation" as given by eqn.(285). One 
then obtains : 



< Ho > = h(M) - 3A/ 2 (M) + i mV + A cr 4 (300) 
We thus derive the LO-effective potential of NGAS as given by 

U (a) = i mV + A cr 4 + h{M) - 3A/ 2 (M), (301) 

where, it is implicitly understood that the "gap-equation" , eqn.(285) is to 
be first solved to obtain M 2 as a function of 'a '. 

One can next carry out the renormalisation programme ( in the LO) by 
noting that 154 : 

(i) the vacuum-configuration corresponds to the absolute (global) minimum 
of Uo(a) , i.e. by solving : 

dU . d 2 U 0l , \ 

d^L = ; ^ U > (302) 

(ii) the renormalised mass in LO is given by : 

2 d 2 U 

m R — k = o-o 5 (ovo) 

(iii) the LO-renormalised coupling strength is likewise defined to be: 

A « = 3T^I— • < 304 » 

where o"o corresponds to the vacuum-configuration as defined by eqn.(302). 
It is directly verified by minimisation of Uq(o~) ( see, eqn.(302)), that the 
global minimum of the former occurs at o"o = , which is consistent with 
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eqn.(288) as it should be. Next, evaluating eqn.(303) at a = one gets 
the renormalised mass : 



m \ = m 2 + l2\I (m 2 R ) = M 2 ( A, a 2 = 0), (305) 

Similarly, after a straight forward calculation, one obtains ref.[154] the 
renormalised coupling as given by: 

*- - 

where again I-i(x) is defined as per the general definition given in 
eqn.(299). 

At this stage several remarks / observations are in order : 

(a) The results contained in eqs. (285,286 ) and eqs. (301-306) were first 
derived by Stevenson 154 in the context of the "Gaussian effective potential 
(GEP)" for the symmetric A0 4 theory and obtained by variational calcula- 
tion using a Gaussian-trial wave-function. 

The reproduction of the results of the GEP in ref. [154] in the LO of NGAS 
demonstrates that the GEP (see, eqn.(301)) is contained in the NGAS as the 
leading order approximation. 

(b) The development of the IPT ( see, chapter-10 ) has therefore, the 
potential to go beyond the Gaussian approximation in a systematic way, or- 
der by order. ( The application of IPT to A0 4 theory falls however, beyond 
the scope of the present dissertation.) 

(c) In view of the equivalence with the GEP in the leading order, all the 
results obtained in the former approximation, are reproduced in the LO of 
the NGAS. In particular, the demonstration of non — triviality of the sym- 
metric \<j) 4 theory in the GEP, being entirely based upon the consequences 
of renormalisation as contained in eqs. (305) and (306), is also reproduced in 
the LO of NGAS. We discuss in the following sub-section, some of the results 
concerning the stability and non-triviality of the theory. 

(d) It must be emphasized, however, that the current scheme, NGAS is based 
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upon entirely different starting assumptions and is much more general 
than the GEP, which is obtained solely due to the choice of the AP as given 
in eqn.(259) and that too, in the leading order. 

We next discuss some of the consequences of the above non-perturbative 
renormalisation scheme obtained in the LO of NGAS leading to the stability 
and non — triviality of the theory. 

11.3 Stability and non-triviality of A0 4 -theory in the LO 
of NGAS 

For the above purpose, it is convenient to start with eqs.(285) and (286) 
, which involve the divergent integral due to the momentum- integration and 
which, therefore, need a suitable method of subtraction. 

Using the subtraction-procedure devised by Stevenson 154 , eqn.(306) can 
be inverted to express A in terms of the observable parameters, \r and 
rriR . This leads to two solutions for A of which, the physical one is given 
by: 

A = [-l/QU(m R )][l + l/2[A il /_ 1 (m i? ) + ] (307) 

( The other solution is A = — (1/2) X R + 0(1/J_i(to.r) . This solution 
can be shown to lead to instability, since the minimum of the EP correspond- 
ing to this solution lies (infinitely) higher than the minimum corresponding 
to eqn.(307)). 

It may be noted at this point that eqn.(307) implies a viable, stable A0 4 
theory when the (unobservable) bare-coupling becomes negative but in- 
finitesimal. This version of the theory has therefore, been designated in 
ref.[155] as precarious A0 4 theory. 

Substituting for A as given by eqn.(307), one can solve for the bare- 
mass, by inverting eqn.(305), after carrying out the subtraction as per the 
Stevenson-prescription. This leads to the following expression for the bare- 
mass : 
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m 2 = m R + 2Io(m R ) / I-i(m R ) + (sub — leading terms) (308) 



With the aid of eqs.(307) and (308), the effective potential, as given 
by eqn.(301) can be recast in manifestly renormalised form involving the 
observable parameters: \ R and m R only. The resulting expression is 
given by : 

U {<J) = U mm + ht m R 2 a 2 - (m R A /l28n 2 )(t - l) 2 - (m R 4 /64n 2 )(t - 1)77, 

(309) 

where 

t = M 2 (a)/m R 2 ; 77 = -An 2 /X R , (310) 

and 

Umin = h{m R ) - ?>\ll(m R ). (311) 
Similarly, the renormalised version of the "gap-equation" is given by : 

(1 -rj)(t- 1) - (16n 2 /m R 2 )a 2 = tint. (312) 

It must be pointed out that, one has to first solve the gap- equation, 
eqn.(312) to obtain t = t(o~) , which is then to be substituted in eqn.(309) 
to infer the cr-dependence of Uq(o~) . 

It may be noted that the gap-equation, eqn.(312) is a transcendental 
equation and its solution exists only when, 

^ 2 < e-Ln = (m R 2 /16n 2 )[e^ + V - 1] (313) 

The domain of validity of the effective potential (EP) is thus restricted by 
the range of a 2 determined by eqn.(313) for any value of n. In particular, 
in the regime of large-coupling (77 — > ), the domain of the EP shrinks 
with i] since cr 2 ^ — > in this limit. This is the situation, therefore, of 
small oscillations about a ~ and it corresponds to the pathological 
situation when \X R \ — > 00. On the other hand, the small coupling regime 
( X R — > 0), which corresponds to n » 1, the domain of EP increases with 
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i] . One can thus summarise that the LO-EP of the symmetric Ac/> 4 theory is 
reasonable and well behaved unless the renormalised coupling is very large. 
One can thus conclude that a non-trivial and stable theory results in the LO 
of NGAS provided the physical coupling Xp is not unusually large. We 
further comment upon the issue of stability of the perturbative vacuum of 
the theory in the following subsection. 

To study the stability issue, it is necessary to compute the effective poten- 
tial based upon the perturbative vacuum (i.e. the vacuum of the free-field 
theory ). This is easily achieved by letting M — > m in all formulae ( see, 
eqn.(264)). Thus starting from eqn. (297) and letting M — > m one obtains, 
after simplification the following expression : 

<H >p = -m 2 a 2 + Act 4 + h + 6Xa 2 I + 3A/ 2 (314) 
2 

In the above, < Tio >p denotes < 0|7io|0 > i.e. VEV of Ti,Q in 
the perturbative vacuum state ; and I n = I n (m 2 ). By definition, the 
effective-potential based upon the perturbative vacuum denoted by U p (a) 
is identified with < Ho > p ■ Hence one obtains : 

U p (a) = -m 2 a 2 + Act 4 + h + 6Xa 2 I + 3A/ 2 (315) 

The renormalised-parameters following from the eqn. (315) are like-wise 
computed and denoted by 

m\ = d 2 U p /da 2 \ a2 = o ; (316) 

Xr = (l/m 4 U P /da%, = o (317) 

where cr 2 = is again the location of the global- minimum of U p (a) , 
as can be readily verified. In this context, it must be emphasized that 
the integrals : I n occuring in eqn. (315) are independent of a . Next, com- 
puting the derivatives of U p (o~) at the minimum, one obtains the following 
expressions for the renormalised parameters based upon the perturbative 
vacuum: 

m\ = m 2 + 12A/ (318) 
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\r = \ (319) 

The requirement of the finiteness of m R and \r then demands that A 
must be — ve(ref.to eqn.(318)), for otherwise fh R would be infinitely large 
since Iq is divergent and the bare ( unobservable ) mass,m 2 > 0. However, 
this would lead to instability since the effective-potential U p (a) will not have 
a lower-bound! This is made manifest by explicitly writing the EP in terms 
of the renormalised parameters: 

U p {a) = l -m\a 2 + \ R a 4 + 3\ R P + h (320) 

To prevent instability of the theory when renormalised about the pertur- 
bative vacuum it, therefore, becomes inescapable that 

\r = A = , (321) 
which is nothing but the triviality scenario ! 

A few remarks/observations are in order, in view of the above results: 

(i) We believe that the result in ( eqn.(321)) constitutes perhaps, the 
most direct demonstration of triviality of symmetric A0 4 theory in physical 
dimensions. 

( ii ) At the same time, the result, eqn.(321) also demonstrates that the 
conclusion of trivialityofthetheoryisanartefactofthenaive perturbation theory 
built and renormalised around the free — field vacuum . As demonstrated 
earlier, the theory renormalised about the NGAS-vacuum leads to a per- 
fectly acceptable, stable and non-trivial A0 4 - theory ( see, eqs. (309-313) and 
discussions following ). 

( iii ) It may be further pointed out that the ground state of the trivial 
theory is still unstable as compared to that in the LO of NGAS, i.e. 

Umin « C- (322) 

which is readily established by referring to eqn.(320) ( with a — ) and 
eqn.(311). 
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This completes our results and discussions regarding the stability and the 
triviality of A0 4 - theory in the context of NGAS in the LO. 



In the next subsection, we discuss the structure of the interacting vac- 
uum in analogy with the results obtained in the case of the AHO/DWO's 
considered earlier. 

11.4 Properties of the Interacting Vacuum State (IVS) 
in NGAS 



The actual/physical vacuum state in presence of interaction is approxi- 
mated in the LO of NGAS by the state: \vac > which is the lowest energy 
state of Ho . In analogy with the results obtained for the AHO/DWO's in 
chapter-9, the structure and properties of this state can be inferred from 
studying the quantum-canonical transformation ( Bogoliubov-Valatin trans- 
formation, ref.[52]) connecting the interacting vacuum state (IVS) with the 
free- field vacuum ( FFV ) state. 

For this purpose it is convenient to start from the Fourier-decomposition 
of the field <f>(x, t) in terms of the free — field creation- and annihilation 
operators analogous to eqn.(266): 



<f>(x,t) = 




where, now 



n k (m) = 2 (2tt) 3 u k (m) ; k° = uj k {m) = ^\k\ 2 + m 2 , (324) 

corresponding to the propagation of the free-field quanta satisfying the 
mass-shell condition : k° 2 — \k\ 2 = m 2 . The free-field operators satisfy 
the standard commutation relations : 

[a(k), a\q)} = Q k {m) 5 3 {k - q). (325) 
Comparision of eqs.(268) and (325) implies that the modified operators : 



B(k)= m 
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A(k) = ^= (326) 

satisfy identical commutation relations : 

[B{k),B\q)\ = 5\k-q) = [A(k), A\q))]. (327) 



It follows, therefore, that the two sets of operators must be connected 
through Bogoliubov-transformation ref.[52] given by : 

B{k) = cosh{a k )A{k) — sinh{a k )A\—k) 

B ] (k) = cosh{a k )A\k) - sinh(a k )A(-k), (328) 

whereas the inverse transformation is given by : 

A(k) = cosh(a k )B(k) + sinh{a k )B\-k) 

A* (k) = cosh(a k )B\k) + sinh{a k )B(-k). (329) 

— * — * 

In the above, a k = f(\k\) , is apriori an arbitrary real function of \k\ , 



i.e. 



a -k = a k = a h ( 33 °) 



However, eqs.( 268 ), (323) and (328) considered together further imply 
that 



exp (2a k ) = — — — = V (331) 
Uk{rn) ^\k\ 2 + m 2 



To show this, consider eqs.(268) and (323) at t — , which can be 
written as: 
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0(x,t) = a + [ f k [ A(k) + A\-k) }e lk " } 



flk(m) 
d 3 k 



/(i k' — » -* ■* 

-==[B{k) + £t(-fc) ]e* ■ * , (332) 
JVt k (M) 



which implies that : 



{ £(£) + B\-k) } 



However, from eqn.(328) it follows that 

{B(jfe) + S^-jfe)} = exp (a fc ){ A(k) + A f (-£) } (334) 
thus leading to the desired result, eqn.(331). 

To obtain the information regarding the particle-content and other fea- 
tures of the IVS it is instructive to first compute the number-density of the 
free- field-quanta residing in the IVS. To this end let us note that the free- 
field-number operator is given by the standard expression : 

N = I ^-^(k)a(k) 

= J d 3 k A\k) A(k) (335) 

Hence the desired number density of the free-field quanta in the IVS is 
given by 

/A , AUk)A(k) . , , 

n(k) = < vac\ — v \vac> , (336) 

where v = spatial- volume of quantisation = / d 3 x. Using eqs.(329), 
eqn.(336) is easily evaluated : 

n$) = ' <B *y ) (337) 
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where a vac (k) is given by eqn.(331) evaluated for M = mn = free- 
particle- mass renormalised about the IVS, \vac >. ( It may be recalled that 
M(a = 0) = rriR and a = define the IVS ). This leads finally to the 
expression : 

/A / 1 w Ukjm) uj k (m R ) 

n{k) = W [ ^y + ^r- 2] - (338) 

To extract further meaningful content from eqn.(338), we note that the bare- 
mass is divergent : (^) ~ 0( ^== = ) where A = momentum cut-off ( see, 
eqn.(309)). Since, according to the standard prescription of the renormalisa- 
tion procedure , the cut-off must be removed ( i.e. A — > oo ) prior to the 
calculation of any physical quantity of the theory , one obtains : 



2 



i 



Sol <^> s **) = < : + ^>~ 5 • < 339 > 

where n(0) = n(/c)| moa; = ( )( ), is the maximum value of n(k) , 

— * 

occuring at = 0. 



Equation (339) provides direct physical content for the non-trivial struc- 
ture of the IVS representing a condensate of off- shell correlated particle-pairs. 
The situation is analogous to the structure of the physical vacuum state in 
case of the hard-sphere Bose-gas 49 and superfluidity 147 . It is therefore , plau- 
sible that eqn.(339) might lead to interesting consequences for T 7^ , as 
happens in the case of the super-fluid and the hard-sphere Bose-gas. 



12. Summary, Conclusions and Outlook 

In summary, a new scheme of approximation in quantum theory, is pre- 
sented which is simple, non-perturbative, self-consistent and systematically 
improvable. The scheme is, in principle, applicable to arbitrary interacting 
systems. We have, however, confined the application of the method to the 
quartic, sextic and octic anharmonic oscillators, to the quartic and sextic 
double well oscillators and to the A0 4 symmetric QFT in the present work. 
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The essential method of this scheme of approximation consists of finding 
a "mapping" which maps the "interacting system" on to an "exactly solv- 
able" model, while preserving the major effects of interaction through the 
self consistency requirement of equal quantum averages of observables in the 
two systems. 

This approximation method has the advantage over the naive perturba- 
tion theory (NPT) and the variational approximation by transcending the 
limitations of both: unlike the variational method, it is systematically im- 
provable through the development of an improved perturbation theory (IPT) 
whereas, in contrast to the case of the NPT, the latter satisfies the necessary 
condition of convergence for all allowed values of the quadratic and the an- 
harmonic coupling strenths ' g ' and 'A'. The method reproduces the results 
obtained by several earlier methods 155 but transcends the limitations of these 
methods in respect of wider applicability, systematic improvement and better 
convergence. 

A remarkable feature of the scheme is that it respects the exact predic- 
tions of super symmetric quantum mechanics (SUSYQM) to a good degree of 
accuracy in case of the sextic anharmonic oscillator and the sextic double well 
potential, when these form a set of "partner potentials". In particular, the 
property of "iso-spectrality" , "positivity"of energy levels and the predictions 
for the "exact" ground state wave function are reproduced with good ac- 
curacy even in the lowest order of approximation. 

We have also investigated the stability properties and the structure 
of th 'effective' vacuum (EV) of the exactly solvable Hamiltonian, H , 
which models the fully interacting system in the leading order. In particular, 
it is shown that the free-field("perturbative") vacuum is unstable for all val- 
ues of the couple strength in comparision with the EV. Moreover, the latter 
is endowed with a rich structure ("dressing") in terms of the free-field quanta 
manifested by the increasing number density of particles with the strength 
of the interaction, in analogy with the case of the super fluid Helium and the 
hard sphere Bose gas. The application of the method to quantum statistics, 
non-oscillator systems and field theory appear to be straight forward. 

The present approach to A0 4 theory reproduces the main results of the 
Gaussian approximation 154 . This is considered quite significant since the 
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two approaches are based upon rather different physical assumptions. How- 
ever, the authors consider the present approach to be more general than 
the Gaussian approximation since the former provides a dynamical expla- 
nation of the latter through the mechanism of altered vacuum-structure in- 
troduced by the interaction. Besides, the authors go beyond the scope of 
the Gaussian-approach in establishing new results, e.g., the calculation of 
the momentum-distribution of the condensate-structure function 'n(k)'. It 
may be emphasized that by going beyond the LO of NGAS, the results of 
the Gaussian approximation can be systematically improved, order-by-order. 

As has been demonstrated, application of NGAS in the LO leads to a 
nontrivial and stable A0 4 theory in the symmetric phase. It is well-known, 
however, that lattice investigations 156 of 4 -field theory indicate the trivial- 
ity scenario and miss the non-trivial version arrived at here and in ref. [154]. 
This result can be succinctly understood as follows: the lattice regularised 
version of the theory corresponds to a finite ultraviolet-cut-off. This means 
that the bare coupling A is small ~ 0(l//^** ce ) and negative for the case 
considered in this work. However, the range of the classical field a remains 
unrestricted. For this reason, there always exist sufficiently large values of 
a (for any given lattice-spacing) such that the term Act 4 occuring in the 
effective potential (see, Eq.(301) dominates over all other terms (Note that 
integrals, I n are all finite on the lattice). Therefore, V latUce (a) becomes 
unbounded from below hence, unstable since 'A' is negative ! 

In a continuum theory, however, the ultra-violet cut-off is never actually 
present: if a cut-off is introduced to regularise the theory, the same has to 
be sent to infinity first prior to considering any other limiting behaviour, 
such as |ct| — > oo . This crucial difference in the order of taking limits: 
(UV-Cut-off) A — > oo and |ct| — > oo, makes all the difference in the physical 
content of the theory in the two approaches and explains why, for any finite 
lattice-spacing, it will not be possible to discover the stable and non-trivial 
version of the theory presented in this work and in ref. [154]. For a detail 
discussion on this important point, see ref. [154]. 

The resulting momentum distribution of the vacuum condensate struc- 
ture function 'n(fc)' deserves special mention as it displays the non-standard 
feature of an appreciable spread in \k\ about the orgin scaled by the renor- 
malised mass of the physical quanta. It is reasonable to expect that this 
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condensate-structure of the physical vacuum persists to finite temperature 
manifesting in observable consequences in the thermodynamic properties of 
the associated system. This would, therefore, constitute a test of the basic 
underlying assumptions of the approximation. 

Finally, it may be worthwhile to compare and comment on related work 
in the recent literarure. As has been remarked earlier, analogous ansatz for 
the ground state and the field-operators derived by Boguliobov transforma- 
tion, has been used in ref. [38] for the case of the anharmonic oscillator. The 
important conclusion that emerges from this study is that a convergent and 
accurate perturbation theory for the energy levels results when the theory is 
developed about the trial vacuum-state. In contrast, the perturbation theory 
is badly divergent 12-14 ' 79 if developed about the non-interacting (perturba- 
tive) vacuum. 

The relation of the present variational approximation (which is equivalent 
in the LO to the Gaussian approximation) to the one-loop approximation 
method 157 has been discussed in detail in ref. [154]. It has been shown that 
the one-loop results are contained, as a special case, in the present results 
(Eqs. (305-306)) in the limit of small bare-coupling, A — > 0+ when 0(A 2 ) 
terms are neglected. This means that the one-loop results remain essen- 
tially perturbative in nature, even though one resums an infinity of ordinary 
Feynmann diagrams at the one-loop order. Besides, the one-loop effective 
potential has a rather restricted domain in both A and a beyond which it 
shows pathological behaviour. In contrast, the effective potential based upon 
the present approximation, even in the LO, has a considerably larger range 
of validity and is genuinely non-perturbative in nature. 

In a spirit similar to the present approach, massless A0 4 theory has been 
variationally investigated in ref. [158]. It has been shown 158 in this work that 
the preferred vacuum is also described by a condensate structure albeit with 
a different momentum distribution function of the condensate particle den- 
sity. In lower (1+1) dimension, the GEP has been derived 55 for the A0 4 
theory employing similar ansatz for the trial vacuum state. However, the 
underlying condensate structure and renormalisation have not been investi- 
gated in these works. 

As remarked earlier, possible applications of the results derived here are 
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envisaged in diverse area of current interest including critical phenomena (in- 
volving a scalar field as the order-parameter 122 ), inflationary cosmology 153 , 
finite temperature field theory 159 , exploration of the vacuum structure 160 of 
pure gluonic-QCD and Higgs sector of the standard model ( by extending 161 
the analysis to the spontaneously broken phase, which corresponds to the 
case of negative bare- mass m 2 < ). 

The present work can be extended in different directions to include : finite 
temperature field theory, quantum-statistics, application to non-oscillator 
systems, super-symmetric theories and quantum field theories involving fermions 
and gauge-fields etc. 
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